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We present the Spitzer/IRS spectra of 157 compact Galactic planetary nebulae 
(PNe). These young PNe provide insight on the effects of dust in early post-AGB 
evolution, before much of the dust is altered or destroyed by the hardening stellar radi- 
ation field. Most of the selected targets have PN-type IRS spectra, while a few turned 
out to be misclassified stars. We inspected the group properties of the PN spectra and 
classified them based on the different dust classes (featureless, or F; carbon-rich dust, 
or CRD; oxygen-rich dust, or ORD; mixed-chemistry dust, or MCD) and subclasses 
(aromatic and aliphatic; crystalline and amorphous). All PNe are characterized by dust 
continuum and more than 80% of the sample shows solid state features above the con- 
tinuum, in contrast with the Magellanic Cloud sample where only ~ 40% of the entire 
sample displays solid state features; this is an indication of the strong link between 
dust properties and metallicity. The Galactic PNe that show solid state features are 
almost equally divided among the CRD, ORD, and MCD. We analyzed dust properties 
together with other PN properties and found that (i) there is an enhancement of MCD 
PNe toward the Galactic center, in agreement with studies of Galactic bulge PNe; (ii) 
CRD PNe could be seen as defining an evolutionary sequence, contrary to the ORD 
and MCD PNe, which are scattered in all evolutionary diagrams; (iii) carbon-rich and 
oxygen-rich grains retain different equilibrium temperatures, as expected from models; 
(iv) ORD PNe are highly asymmetric, i. e. bipolar or bipolar-core, and CRD PNe 
highly symmetric, i.e., round or elliptical; point-symmetry is statistically more common 
in MCD than in other dust class PNe. By comparing the sample of this paper to that 
of Magellanic Cloud PNe we find that the latter sample does not include MCD PNe, 
and the other dust classes are differently populated, with continuity of the fraction of F, 
CRD, ORD, and MCD population from high to low metallicity environments. We also 
find similar sequences for CRD PNe in the Galactic disk and the Magellanic Clouds, 
except that the Magellanic Cloud PNe seem to attain higher dust temperatures at sim- 
ilar evolutionary stages, in agreement with the observational findings of smaller dust 
grains (i.e, lower radiation efficiency) in low metallicity interstellar environments. 

Subject headings: planetary nebulae: general 



1. Introduction 

Planetary nebulae (PNe) are the gas and dust envelopes ejected toward the end of the evolution 
of low- and intermediate-mass stars (LIMS, 1-8 M Q ), at the tip of the thermally-pulsing asymptotic 
giant branch (TP-AGB; e. g., Herwig 2005) phase. The stellar ejecta carries both the products of 



1 based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under a contract with NASA. 
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nucleosynthesis, such as carbon and nitrogen, and the a-elements, such as oxygen, neon, argon, 
and sulfur, whose net yields tend to be near zero in this mass range. 

Dust may be fundamental to PN formation, since radiation pressure on the dust grains formed 
at the surface of AGB stars may be important in triggering the envelope ejection. Planetary nebulae 
are ideal dust probes during their early evolutionary stages, since sputtering of the dust grains is 
expected to affect their dusty nature in their lifetime. Dust particles in PNe are typically cool (50- 
150 K) and radiate in the mid infrared, producing a near-blackbody continuum spectrum (Cohen 
& Barlow 1974) that contributes to ~ 40% of the total emergent flux (Zhang & Kwok 1991), and 
that peaks between 25 and 60 //m. Spitzer IRS spectroscopy has been very successful in tracing 
this dust component in PNe, as well as in a variety of other stellar targets including the AGB 
stars. Previous systematic PN studies based on IRS spectra include Galactic bulge (GO program 
3633, Bobrowsky), halo (GO program 20049, Kwitter), disk (GTO program 40035, Bernard-Salas), 
and Magellanic Cloud (GTO program 103, Houck, and GO program 20443, Stanghellini) PNe. 
These data collectively show PN spectra with thermal dust continua, nebular emission lines, and 
a variety of dust signatures, with characteristics of carbon-rich and oxygen-rich compounds and 
different types of grains. From the comparative study of the different samples some trends have 
emerged: it appears that oxygen dust features, such as crystalline silicates, are more common in 
the bulge population (Gutenkunst et al. 2008; Perea-Calderon et al. 2009) than in the Galactic 
or Magellanic Cloud PNe (Stanghellini et al. 2007, hereafter S07; Bernard-Salas et al. 2009). In 
particular, in the Magellanic Cloud PNe the oxygen and carbon dust features are never observed 
in the same spectrum. It is apparent that the galaxy metallicity has a strong effect on the nature 
of dust (S07; Bernard-Salas et al. 2009), and, from the Magellanic Cloud sample, it is apparent 
that planetary nebula bipolarity correlates strongly with oxygen-rich dust PNe (S07). Finally, gas 
and dust chemistry in Magellanic Cloud PNe are strongly correlated (S07). 

What was lacking from the Spitzer database was a large, homogeneous sample of Galactic PN 
spectra, especially for young, unevolved PNe, to study the early dust features and their evolution 
across the Galactic disk. In this paper we present the study of 157 Galactic angularly small targets 
classified as PNe, observed during the last cryogenic cycle of Spitzer. We collected the present 
sample to gain insight on the impact of dust and metallicity in stellar evolution, a fundamental 
question in astrophysics. The extensive target list is essential to explore the different dust types with 
statistical significance, and to study them across the Galactic disk, and in relation to the Galactic 
populations. The large sample presented here fills an important gap in the Spitzer program. Our 
compact Galactic PN dataset is also complemented by a WFC3/ HST imaging survey of about 
one-third of the targets presented here. The morphological analysis, and that of the central stars 
(CS) from the WFC3 data will complement the dust analysis presented here, allowing a much more 
detailed view on how dust affects the post- AGB and PN evolution in our galaxy compared to other 
environments. 

Systematic comparison of IRS spectra of PNe belonging to different populations such as the 
Galactic disk, halo, bulge, and the satellite galaxies, can address many key questions that have 
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opened up in recent years regarding AGB stars as well. First, the absence of heavily obscured AGB 
stars in the Magellanic Clouds (Trams et al. 1999; Groenewegen et al. 2000; Garcfa-Hernandez et 
al. 2009), in contrast with their Galactic counterparts (e. g., Garcfa-Hernandez et al. 2006, 2007a), 
seems to indicate that, on average, lower metallicity environments such as those of the Magellanic 
Clouds are less favorable to dust production. The relative number of C-rich versus O-rich AGB 
stars in galaxies of the local group increases with decreasing metallicity (Cioni & Habing 2003; Cioni 
et al. 2003; Schultheis et al. 2004), which shows that mass-loss efficiency depends on metallicity 
(Hofner 2011). 

In this paper we present the 5-40 /um IRS spectra from Spitzer program GO 50261. In §2 
we describe the target selection, exposures, and observing strategy. The data analysis, including 
spectral extraction, classification of the dust types, and continuum fitting, are presented in §3. In 
§4 we explore the relation of the dust characteristics with respect to other physical properties of 
the PNe in our sample, and in §5 we extend the comparison to PNe in the Magellanic Clouds. The 
summary and conclusions are in §6. 



2. Observations 

The aim of the target selection is to cover a homogenous and as much as possible complete set 
of compact Galactic PNe, most of which should be young. To this end, we selected PNe smaller 
than 4" in apparent size. Galactic PNe with diameters < 4" should be younger than ~ 2.0 x 10 3 
yr (Villaver et al. 2002a, b) if they have heliocentric distances smaller than 6 kpc and if their 
expansion velocities are typical (20-40 km s" 1 ). Statistical distances (Stanghellini et al. 2008) 
help us distinguish which PNe are nearby and young from those that are distant and evolved. 
The population of spectroscopically confirmed Galactic PNe is listed either in the Strasbourg-ESO 
catalog of Galactic PNe (Acker et al. 1992) or in the MASH survey (Parker et al. 2006). Of the 1143 
PN in the Strasbourg-ESO catalog, 143 are point sources, and 86 have 6 < 4". The MASH survey 
gives another 2 PN with 9 < 4" (but no point sources). From the 230 PNe thus selected we explicitly 
eliminated the PNe already observed with Spitzer/IRS. Most of the exclusions are Galactic bulge 
and halo PNe, only 10 compact PNe in the Galactic disk had been observed with Spitzer before 
with similar observing configuration. Finally, we also exclude a dozen targets previously observed 
with IRAS or MSX. This filtering yielded a final sample of 157 compact Galactic PNe, representing 
for the most part the disk population. 

Twenty-six PNe in our sample may actually be bulge PNe. The fact that a PN belongs to 
the bulge rather than the disk is controversial, if nothing else for the reason that the Galactic 
PN distances are not known with high precision, thus finding a PN in the direction of the bulge 
might not be sufficient to exclude it form the disk population. Stanghellini & Haywood (2010) 
made this selection for all Galactic PNe based on the best distance scale available, and included 
a brightness criterion, but uncertainties in the separation of the bulge to disk PN samples are of 
course possible. Furthermore, many of the PNe of our sample do not have a distance determined 
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Fig. 1. — The spatial distribution of our targets (filled circles) compared to the general PN popula- 
tion in the Galaxy (open circles). This plot is limited to those PNe whose statistical distance could 
be calculated, thus it represents subsamples of both our target sample and the general Galactic PN 
population. 
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at all, given the lack of measurable angular diameter. Knowing the radial velocities of the targets 
could help, unfortunately at this time velocities are available only for a minority of the sample, and 
none for the PNe whose distances are not known. 

Three PNe in our sample could belong to the Galactic halo. In Table 1 we list the observed 
targets with their IAU designation (column 1), their common name (column 2), and their equatorial 
coordinates (equinox 2000, columns 3 and 4). In columns 5 through 7 of Table 1 we list the IRS 
campaign number, the observing mode, and the exposure time respectively. In the observing mode 
column, "0" is for the combination of SL (short-low, 5-14 jum) and LL (long-low, 14-40 /urn) modules, 
and "1" for the combination of the SL, SH (short-high, 10-20 fim) and LH (long-high, 20-40 fxm). 
Nearly all targets were observed with the SL module and "0" and "1" denote respectively whether 
the longer wavelength part of the spectra where observed at low or high resolution. 

In Figure 1 we show the space distribution of our targets whose distance is known (92 PNe) 
compared with the general PN population (Stanghellini &; Haywood 2010), where their distance 
from the Galactic plane is plotted against their radial distance from the Galactic center, and where 
we assume that the Galactic center is 8 kpc from the sun. It is clear from this plot that the selected 
sample well represents the Galactic disk population of PNe. Most of the PNe in this study are 
actually young, being selected to be compact in apparent size. 

We acquired, reduced, and analyzed the Galactic PN IRS/Spitzer spectra in a similar manner 
to that for Magellanic Cloud PNe (S07). Our AOR were built with the following requirements: 
(1) we aimed at obtaining 5-40 ^m spectra, with different modes depending on the estimated PN 
brightness; (2) we estimated the fluxes of our targets based on the upper limit provided by IRAS 
(Zhang &; Kwok 1991). We chose the observing mode, computed the exposure time and the number 
of cycles needed for each mode using Spitzer IRS PET, the saturation levels for all modes, and 
a conservative choice for the fluxes. For the PN with very low fluxes, based on the IRAS upper 
limits and our experience with IRS spectra of Magellanic Cloud PN, we observed in SL and LL 
modes, with 3 integrations of respectively 14 and 30 seconds each. Peak up observations were not 
necessary, since we derive all coordinates form the astrometric catalog by Kerber et al. (2003), 
where the precision is better than 0.35" . 

3. Data analysis: calibration, spectral extraction, and dust classification. 

Raw data were retrieved from the Spitzer archive together with the latest pipeline calibration 
files. Initial file preparation included removing Bowing of spectral orders in the LH data using 
DARK SETTLE; cleaning rogue pixels using IRSCLEAN MASK, where we used custom masks for 
rogue pixels cleaning for all our data; and subtracting the two nod positions for the SL and LL 
data to obtain sky-subtracted spectra. All low- resolution spectral images were checked free from 
background source contamination. We did not subtract the sky in the high resolution data, but 
the sky background is very low for these bright sources. We then extracted the ID spectra from 
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Fig. 2. — The spectral templates of dust types found in our compact Galactic PN sample. 
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the co-added 2D spectral images using SPICE version 2.2.. 

The spectral merging, averaging, and continuum fitting was done with the package SMART 
(Higdon et al. 2004). The final spectrum for each dataset includes merging, cleaning and sub- 
tracting any spurious jump, and averaging the nod positions and orders. We subtracted a spurious 
pedestal above 14 fim from the spectrum of PN G012.5-09.8, so both final reduced spectrum and 
the fitted black body curve have been modified accordingly. We also subtracted a pedestal along 
the whole spectrum of PN G041. 8+04.4 and PN G044. 1+05.8 to bring it down to a level of 0.0 Jy 
at 5 //m. 

Three distinct components are apparent in the spectral energy distributions: the dust con- 
tinuum, the nebular emission in the form of collisionally excited lines of atomic gas and, in the 
majority of the targets, solid-state dust emission features. An analysis of the abundances in the 
nebular atomic gas will be presented in a future paper; for now we focus on the spectral signatures 
of the dust. We fit the dust continuum in each target with the black-body fitting routine described 
in the package SMART/IDEA in order to determine a characteristic dust temperature and IR lu- 
minosity. We performed a continuum fit on each spectrum after masking the nebular and solid 
state emission features and the low-signal region long- ward of ~ 37^m where the system sensitivity 
drops rapidly. The temperature and luminosity were varied to minimize the RMS deviation, and 
the IR luminosity was derived from the integral of the Planck function at that temperature over 
all wavelengths. In most cases the continuum was well sampled over the bulk of the observed 
spectrum, through the whole SED, and we constrained the fitted function to pass through or below 
the average continuum. 

The form of the fitting included an emissivity term, r oc A _a , (a is given in column (5) of table 3 
for the converging fits) and the fitting routine finds a solution for both alpha and the proportionality 
constant (see also Hony et al. 2002). The black body fits were robustly obtained for about half 
of the PN sample, and in these cases the formal uncertainties of the fit to temperature and IR 
luminosity were small, owing to the very slow change in the Planck function with wavelength. For 
cases where the fits did not converge with the proposed method (flagged as N in Table 3, column 4), 
we did not use the derived dust temperatures in the tables or Figures of this paper. We plan to use 
more sophisticated fitting models in a future analysis. For the PNe whose black body fit converge, 
we have validated them further by comparing the fits against the the 65 //m fluxes from the Akari 
data archive, and the 60 //m fluxes from IRAS archivqj in order to determine whether the IRS 
spectra do sample their flux maxima (i.e., the 60 [or 65] /im flux is lower than the maximum flux of 
the IRS spectra). In these cases we also have checked whether the 60-65 p fluxes in the literature 



2 In Table 3, column (6) we give the satellite fluxes we have used to asses the validity of the fits, where we used 
the Akari 65 /im flux (Yamamura et al., 20f0) when available, otherwise the 60 /im IRAS flux (given in Acker et al. 
f992). Errors are given when available from the Akari data archives; flux errors are not available for IRAS, only flux 
quality, where a low quality indicates an upper limit to the 60 fim flux, which is still valid, in our case, to constraint 
the continuum fit. 
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were not lower that our continuum fit extrapolations at the same wavelength. The majority of the 
(converging) fitted continua are perfectly compatible with the long wavelength photometry where 
available (55 PNe, flagged with A in Table 3, column 4); the IR flux long- ward of 60 /im shows that 
the observed SEDs sample the functional maximum of the fits. Furthermore, our extrapolated fluxes 
and those in the literature correlate very well with one another, with a global (linear) correlation 
coefficient of 0.82. We use the A fits with confidence in our group analysis and figures. 

Another 14 PNe are below the limit of detectability of the Akari and IRAS satellites at 60-65 
jum, although their continuum fits seem to be excellent; they are flagged with B in column (4) of 
Table 3. 

Finally, 10 PNe have (converging) black body fits which are badly constrained by the 60-65 //m 
archival flux (they are flagged with C in Table 3, column 4). What happens in these cases is that 
the 60 and/or 65 ^m fluxes are considerably higher than the IRS flux maxima. This discrepancy 
can have several causes, including (a) a bad fit of the continuum, (2) the presence of other sources, 
or (3) high background at 60-65 fim. We have inspected these sources and found 4 possible bulge 
PNe among them, which could imply higher than average background, and one PN (PN G356.5- 
03.6) whose WISE data archive inspection within the Spitzer beam disclosed multiple sources which 
could cause the mismatch between the extrapolated fits and the 60-65 /xm archival fluxes. For the 
remaining PNe there is not a good explanation for the flux mismatch, likely an unexplained high 
background at long wavelength, and they will be studied in more detail in the future. To be very 
conservative with our conclusions we will not use the black body fits that do not comply with the 
available 60-65 ^m photometry in the Figures (29-36) or analysis (Section 4) of this paper. 

While these black-body fits reflect the physical origin of the continuum shape, this character- 
ization is only broadly representative of the thermal dust. Indeed, the dust within a given nebula 
may be composed of a variety of chemical species, each with their own emissivity, at varying dis- 
tances from the source of excitation. Thus, there may in reality be an ensemble of thermal emission 
within the nebula which is not perfectly represented by a single black-body temperature. We revisit 
this point and the impact it could have in the analysis presented in Section 4.2. In addition, the 
residuals in the black-body subtracted continuum may complicate the interpretation of very broad, 
weak, solid-state features. For those, one needs to perform polynomial fits and subtractions (see, 
for example, Garcfa-Hernandez et al. 2010). 

In Table 2 we give a description of our dust classification scheme. We define therein the 
dust classes, and also several subclasses that we encounter in the analysis of our large sample of 
compact Galactic PNe. As for the Magellanic Cloud sample, we classify the spectra based on 
their dust features. A detailed analysis of the atomic emission lines and elemental abundances 
will be published elsewhere. We distinguish among (1) the absence of prominent dust features, 
(2) the presence of carbon-rich dust features, and (3) the presence of oxygen-rich dust features. 
PNe lacking prominent molecular /dust emission features and with very low continuum are defined 
as F. Those with carbon-rich features are in the CRD class, with subclasses determining weather 
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Fig. 3. — Featureless spectra, in PN G order. Dotted lines correspond to black body continuum 
fits of type A or B (see text, and Table 3). 




Fig. 4. — Featureless spectra, cont. 




Fig. 6. — Featureless spectra, cont. 
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they show aromatic and/or aliphatic features. Similarly for the oxygen-rich dust features, or ORD, 
where the subclasses indicate the presence of crystalline and/or amorphous dust grains. Finally in 
many cases the spectra show both CRD and ORD characteristics, and we call these dual-chemistry 
PNe as MCD (or mixed-chemistry dust) PNe. In Figure 2 we present the array of template spectra 
according to dust class and subclass. Table 3 lists, for all targets, the dust class (column 2), and 
subclass (column 3). 

In Figures 3 through 26 we plot all IRS PN spectra from our program, and we overplot (as 
dotted lines) the continuum fits for A and B fit types (see Table 3). There are a total of 25 PNe 
with continuum-only (F) spectra; one is shown in Figure 2 (top panel) and the remaining 24 in 
Figures 3 through 6. None of these spectra shows very obvious molecular /dust emission features 
above the continuum, and the continua are very faint. In a couple of cases (PN G038.4-0.3. and 
PN G038.7-03.3, both in Fig. 4) the dust profile has a different shape than the other F PNe, in 
particular there is a suggestion of a bump in the ~20-30 /um spectral region, thus the F designation 
might be uncertain. Actually, there may be some ORD or CRD PNe among the subsample of 
objects with featureless spectra but the weakness of the dust features prevents to confirm the 
dominant chemistry (ORD or CRD) of their circumstellar shells. It is worth noting that two of the 
F spectra (PN G356.5-03.6 and PN G359.3+03.6, both in Fig. 6) show a jump between the SH and 
LH modules (~ A < 19.5/xm) that, in principle, may be attributed to an extended nature of these 
sources at the longer wavelengths. However, due to the compact nature (<4") of our PN sample 
and because we do not see any mismatch between the SL and SH modules, we suspect that this 
LH flux excess is possibly due to a known systematic effect that may affects some LH observations 
(see e.g., Garcia-Hernandez et al. 2007b) or, less likely, to background emission. In general the 
background is very faint relative to the source in these LH observations, but there is really no 
way to check if this is indeed the case. In any case, this slight mismatch between the SH and LH 
observations does not affect any of our results. In these case, we have scaled the continuum for the 
final fits. 

Carbon-rich dust features can be aromatic or aliphatic. The aromatic features are represented 
by the well-known family of infrared features at 6.2, 7.7, 8.6, and 11.2 fim usually attributed to 
polycyclic aromatic hydrocarbons (PAHs; e.g., Leger & Puget 1984), or, as more appropriately 
defined by Kwok & Zhang (2011), to aromatic infrared bands (heretofore, AIBs). The AIBs have 
been observed in several types of astronomical sources. It is worth noting that in several cases a 
single (and usually very weak) 11.2 /im feature in the spectra, especially if the feature is narrower 
than common AIB features, is due to interstellar rather than circumstellar emission, thus the nebula 
would not be classified as CRD solely on this emission feature. In Figure 2 we show a typical 
example of obvious AIB features in PN G068.7+14.8, where all 4 bands are clearly present. PNe 
with spectra showing any of these bands are classified as CRD; our sample includes 12 additional 
such spectra. 

A larger number of the CRD spectra show other emission features that are likely related 
to carbon-rich dust as well, but of different types. The 6-9 /im feature can be attributed to 
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Fig. 7. — CRD spectra, in PN G order. Both aromatic, aliphatic, and aromatic/aliphatic spectra 
are shown here, see Table 3 for complete list. On the left panels we show the complete spectra, 
while right panels show the 5-12 fim sections of the spectra, where the 6.2, 7.7, 8.6, and 11.2 ^m 
AIB band positions have been indicated with vertical lines. Dotted lines have the same meaning 
than in Fig. 3 



-16- 



0.8 — ' — ' — ' — 1 — ' — I — ■ — 1 — ■ — ■ — ■ — 1 — ' — 1 — 1 — 1 — 0.5 f 




10 20 30 40 5 10 15 

X [firr\] X [firry] 



Fig. 8. — CRD spectra, cont. 
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Fig. 9. — CRD spectra, cont. 
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Fig. 10. — CRD spectra, cont. 
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Fig. 11. — CRD spectra, cont. 
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hydrogenated amorphous carbons (HACs), large PAH clusters, or very small carbon grains. On the 
other hand, the broad 10-15 /im emission (centered at about 11.5 /im) and the 25-35 /im emission 
(the so-called 30/iin feature) are usually attributed to SiC (e.g., Speck et al. 2009) and MgS (e.g., 
Hony et al. 2002), respectively. However, the observed variation of these broad features is quite 
consistent with the variety of properties of HACs, which are manifest in a wide range of different 
spectra depending on their physical and chemical properties (e.g., size, shape, hydrogenation; Scott 
& Duley 1996; Scott et al. 1997; Grishko et al. 2001). In particular, Grishko et al. 2001 showed 
that HACs can explain the 21, 26, and 30 /im features. Zhang et al. (2009) also argue that MgS 
is very unlikely the carrier of the aliphatic 30 /im emission seen in C-rich evolved stars. Thus, it 
is very unlikely that the carriers of the broad 11.5 um and 30 /im are related with SiC and MgS, 
respectively, as it was suggested in the past. We tentatively define all these features as aliphatic 
dust features, found in 22 PNe of our sample. 

Three of the CRD spectra in our sample show both aromatic and aliphatic features, thus the 
total number of CRD PNe in our sample amounts to 38, or ^25% of the PN spectra. All CRD 
spectra are shown in Figures 7 through 12, where we show complete spectra in the left panels of 
each figure, and the 5-15 /im section of the spectra in the right panels, to characterize the presence 
(or absence) of AIB features. Most of the CRD spectra are predominantly aromatic or aliphatic, as 
noted in column (3) of Table 3. In a few cases both sets of features are present, and we mark the 
PN as aromatic/aliphatic type. In a minority of cases the CRD dust features are very faint, and 
had the spectra been shallower we would have classified them as F. Their dust type in column (3) 
of Table 3 is then flagged as uncertain. Interestingly, two CRD aliphatic PNe, PN G063.8-03.3 and 
PN G235.3-03.9, display strong Cqq fullerene features (Garcfa-Hernandez et al. 2010), and are part 
of a small group of Galactic and Magellanic Cloud PNe that show the infrared emission bands of 
this complex molecule (Cami et al. 2010; Garcfa-Hernandez et al. 2010, 2011b). PNe that contain 
fullerenes are typically low-excitation and low-mass C-rich objects displaying the broad 11.5 and 
30 /im aliphatic features. The detection of fullerenes in the H-rich circumstellar shells of these PNe 
(see also Garcfa-Hernandez et al. 2011a) has challenged our understanding of the formation of 
these complex fullerene molecules in space, indicating that fullerenes may be more abundant than 
previously believed and demonstrating that normal C-rich PNe are important factories of complex 
organic compounds. 

Carbon-rich dust features have been found in 42 additional PNe, where crystalline silicate 
bands (e.g., Waters et al. 1998) are also present. These objects, mixed-chemistry dust (or MCD) 
PNe, showing both AIBs (C-rich) and crystalline silicates (O-rich) in their shells (e.g., Perea- 
Calderon et al. 2009), have not been found in the Magellanic Clouds to date (S07; Bernard-Salas 
et al. 2009), but are very common in the Galactic Bulge (Gutenkunst et al. 2008; Perea-Calderon 
et al. 2009). A template spectrum of MCD type PNe is shown in Figure 2, while all other MCD PN 
spectra are presented in Figures 13 to 19, where we show complete spectra in left panels, and the 
5 < A < 15/im spectra in the right panels, as in the CRD PN spectra. In the MCD PN plots we also 
indicate the prominent crystalline silicate features at 23.5, 27.5, and 33.8 /im - usually attributed to 
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olivine and pyroxenes - with vertical dotted lines. Crystalline silicates and other oxygen-rich dust 
features (e.g., the amorphous silicate bands at ~10 and 18 /um), not associated with carbon-rich 
dust, have been detected in 45 PNe. Of these, 16 have only crystalline silicate features, 24 show 
possible amorphous oxygen features, and five show possible signatures of both types (amorphous 
and crystalline) of oxygen-rich dust. We call all these objects ORD (oxygen-rich dust) PNe. The 
details of their subtypes are given in Table 2, and their spectra are shown in Figures 20 to 26. 

It turns out that 7 of the 157 targets acquired by us are misclassified as PNe, thus their spectral 
analysis will be published elsewhere. The IRS spectra of the remaining 150 Galactic PNe indicate 
that many more show dust features than their Magellanic Cloud counterparts, as expected, given 
their relative metallicities (S07). Only ~17% of the Galactic sample does not show molecular /dust 
features, while the featureless PNe in the Clouds represented >41% of the samples (S07; Bernard- 
Salas et al. 2009). The PNe spectra with molecular /dust features in the Galaxy are distributed 
evenly among the ORD, CRD, and MCD types, respectively representing approximately 30, 25, and 
28 % of the sample; in the Magellanic Clouds, on the other hand, most of the dust-rich PN spectra 
are CRD, and only very few PNe show ORD features, while MCD PNe have not been detected 
in the Magellanic Clouds. A more detailed analysis (e.g., integrated fluxes, positions, etc.) of the 
organic (AIBs) and solid-state features (aliphatic dust, amorphous and crystalline silicates) seen in 
our Galactic Disk PN sample in comparison those of the Magellanic Clouds will be presented in a 
forthcoming paper (Garcia-Hernandez et al. in preparation). 

4. Correlations between IR spectral characteristics and other nebular and stellar 

properties 

4.1. Ancillary data from the literature 

From the IRS analysis in this paper we have acquired the knowledge of several IR-related 
properties of the PNe. In Table 3, together with the class and subclass of the dust spectra, we 
give the IRS-spectra related parameters, namely, the fit type (column 4), the emissivity exponent 
(column 5), the 60-65 /xm flux from the Akari or IRAS surveys, used to assess the validity of 
the continuum fit (column 6) the dust temperature and its uncertainty (column 7), the total IR 
luminosity (column 8), and the infrared excess, derived classically as in Pottasch (1984) (column 
9). Dust temperatures are given only in the cases where the black body fit converge, either with 
simple black-body or with an emissivity term. It is worth noting that the different dust types 
show a difference in emissivity, as expected, with the average of the exponent a in CRD PNe being 
0.52±0.73, compared to an average value of 3.16±2.06 in ORD PNe. The infrared luminosity was 
determined by integrating the flux under the black body fit of the IRS dust continua, then scaled 
to absolute luminosity by using the statistical distances. The PNe whose fits converge may have 
A, B or C fit type (column 4), as described in section 3. In Figures 29 and beyond we plot Td us t 
and Lir only in the case of an A or B black body fit type. 
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Fig. 13. — MCD spectra, in PN G order. On the left panels we show the complete spectra, where 
the crystalline silicate bands at 23.5, 27.5, and 33.8 /xm have been marked with vertical dotted 
lines; right panels show the 5-12 /im sections of the spectra, where the 6.2, 7.7, 8.6, and 11.2 fj,m 
AIB band positions have been indicated with vertical lines. Dotted lines have the same meaning 
as in Fig. 3 
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Fig. 14. — MCD spectra, cont. 




Fig. 15. — MCD spectra, cont. 
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Fig. 16. — MCD spectra, cont. 
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Fig. 17.— MCD spectra, cont. 
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Fig. 19. — MCD spectra, cont. 
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Other properties of the PNe that we have used in this study have been listed in Table 4. 
The most important parameter to gain astrophysical insight from the analysis of PN properties 
is of course their distances. Individual heliocentric distances determined with reliable methods, 
such as cluster membership, are not available for most of the PNe in our sample. We then use 
the best statistical distance scale available, based on the Stanghellini et al. (2008) Magellanic 
Cloud calibration. In order to get statistical distances on this calibrating scale we need the angular 
diameters, and the 5 GHz (or H/3) fluxes. The fluxes, and the upper limits to angular diameters, are 
available for all targets, but actual angular diameters are known only for 50 of the 150 confirmed 
PNe in our IRS sample, thus their distances are readily available only for approximately a third 
of the sample (see also Stanghellini & Haywood, 2010). An additional 42 PNe have been observed 
by us with WFC3 (Shaw et al. in preparation), and from these space-based images we were able 
to derive the angular sizes, thus calculate their distances with the statistical scheme described 
above. Distances from the Galactic center, and the PN sizes derived as described in Stanghellini &; 
Haywood (2010), are listed in Table 4. 

Plasma diagnostics for the observed PN based on optical spectra are scarce in the literature, 
with only a dozen PNe having adequate electron density and elemental abundance. We have used 
the [SII] A6717-6731 emission lines measured from low resolution spectroscopy and published by 
Acker et al. (1992) to derive an estimate of the nebular densities, by assuming electron temperatures 
of 10,000 K and by using the nebular routines in IRAF/stsdas (Shaw & Dufour 1995). The densities 
calculated are listed in column (4) of Table 4. There is not enough information to determine oxygen 
or other abundances for most of the compact PNe in our sample. The abundance analysis based 
on the IR lines in the IRS spectra will be the subject of a future paper. 

For all other nebular parameters, we used the Stanghellini & Haywood (2010) recent, complete, 
and uniform collection of Galactic PN data. 



4.2. Segregation and evolution of dust-types 

In Figure 27 we show the distribution of dust types as projected on the Galactic plane, where 
Xq=D cos(b) cos(l), Yq=D cos(b) sin(l), D is the heliocentric distance, 1 and b are respectively the 
Galactic longitude and latitude. Diamonds represent the CRD, squares the ORD, triangles the F, 
and circles the MCD PNe. It appears that there is mild segregation of dust type within the Galactic 
distribution, in particular the MCD PNe tend to be concentrated within the plot. To explore this 
segregation further in Figure 28 we show the histograms of dust types as they are distributed in 
relation of their distance from the Galactic center, and we recover the accumulation of MCD PNe 
with Rq < 9 kpc, unlike the other dust types. CRD and ORD PNe seem to be rather uniformly 
distributed. 

Of the three (possible) halo PNe in our sample, two belong to the CRD class, both with 
aliphatic features, while one shows ORD features. We do not observe a spatial segregation of the 
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Fig. 20. — ORD spectra, in PN G order. On the left panels we show the complete spectra, where 
the crystalline silicate bands at 23.5, 27.5, and 33.8 fim have been marked with vertical dotted 
lines; right panels show the 5-12 fim sections of the spectra. Dotted lines as in Fig. 3 
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Fig. 21. — ORD spectra, cont. 
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Fig. 22. — ORD spectra, cont. 
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Fig. 24. — ORD spectra, cont. 
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Fig. 25. — ORD spectra, cont. 
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Fig. 26. — ORD spectra, cont. 
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CRD the dust subtypes. On the other hand, by examining the latitude distribution of the ORD 
subtypes we find that the crystalline ORD PNe have typically lower Galactic latitude (< |b| >= 
2.95 ± 1.76) than the amorphous ORD PNe (< |b| >= 6.33 ± 3.51). This is consistent with the 
former type of ORD PNe being the progeny of high mass AGB stars, while the latter are possibly 
the remnants of the lowest mass progenitors. 

It is worth noting that PN samples in the Galactic bulge, such as those studied by Gutenkunst 
et al. (2008) and Perea-Calderon et al. (2009) show a majority of dual chemistry PNe, unlike 
other known samples. Our Galactic disk sample broadly selects against bulge PNe, but the mere 
definition of Galactic Bulge PN (e.g., Stanghellini & Haywood 2010) is affected by the PN apparent 
angular radii, which can be very uncertain for compact PNe. By using the definition of bulge PN 
by Stanghellini & Haywood (2010) we determine that 26 of our 150 PNe might belong to the bulge. 
We will discuss in more detail the dust properties of bulge, disk, and Magellanic Cloud PNe. In the 
following analysis, and in all subsequent plots, we exclude potential bulge PNe from the Galactic 
sample. 

In Figure 29 we show the dust temperature plotted against the linear size of the PNe in our 
sample whose dust temperature and distances are available. Different type PNe are plotted with 
different symbols (see Fig. 27 legend), and temperature error bars are also plotted, while a typical 
error bar for the radii is given in the plot. At zeroth order, if all PN progenitors were of identical 
mass and metallicity, and if the PNe expand at constant expansion velocity, one would expect log 
Tdust -0-4 log Rpn if the dust grains do not evolve in size (see e. g. Lenzuni et al. 1989). In 
the snapshot of evolutionary stages provided by a sizable observational sample we see the scatter 
due to the different dust types, masses, metallicity, shell acceleration, and alternative evolutionary 
paths. It appears that the CRD PNe are concentrated toward the high dust tempeatures across a 
wide domain of radii. Furthermore, there is a rather tight correlation between dust temperatures 
and physical radii of CRD PNe, with correlation coefficient of -0.99, which is suggestive of a rather 
narrow initial mass and metallicity range. The slope of the CRD PNe on this plane is shallower 
than that found for non-evolving grains (Lenzuni et al. 1989), suggesting there is some type of 
evolution of the dust grains in CRD PNe. 

The solid line in Fig. 29 represents the least squares fit for the CRD PNe, log Td us t=-0.250 log 
Rpn+6.41, where the temperatures are measured in K and the nebular sizes in cm. The ORD PNe 
do not follow a narrow sequence in this plane. A possible explanation for the markedly different 
domains of dust temperatures in CRD and the bulk of ORD and MCD PNe could reside in the 
different heat capacity of dust grains as function of grain size and composition. Li & Draine 
(2001) have shown that equilibrium temperatures of small carbonaceous grains are generally higher 
than for silicate grains in the ISM. PN dust should behave similarly to the ISM one, as the ISM is 
mostly composed by recycled dust from the evolution of LIMS, thus this could explain the observed 
distribution of temperatures of the CRD vs. the other dust type PNe. 

In Figure 30 we show the dust temperature versus the IR luminosity (in solar units) as deter- 
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Fig. 27. — Galactic X and Y location for F (triangles), CRD (diamonds), ORD (squares), and 
MCD (circles) PNe. It appears that most MCD PNe are more centrally concentrated than other 
types 
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Fig. 28. — Distribution of the distances to the galactic center for CRD (top), ORD (middle), and 
MCD (bottom) dust PNe. 
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Fig. 29. — Dust temperature plotted against linear nebular radii (in cm); Galactic disk sample. 
Symbols are like in Fig. 27. The solid line corresponds to the best fit to the CRD PNe. The 
error bars for T— dust have been derived from the continuum fits. The error bar on the left corner 
indicates a reasonable guess derived from the statistical distance uncertainty, corresponding to 
Ad=0.26. 
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Fig. 30. — Dust temperature vs. the integrated IR luminosity. Symbols are as in Figure 27. The 
error bars have the same meanings as in Fig. 29. 
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mined from the black body fits of the continua. In this plot the CRD PNe are found in a sequence, 
where PNe characterized by hot dust are more luminous in the IR. On the other hand, ORD and 
MCD PNe are found scattered on this plane. If for example we look at PNe with log Lir/L q >3 
we find several very hot CRD PNe, while both ORD and MCD PNe have low/intermediate dust 
temperatures. This effect could be dominated by the different heat capacity of the carbon-rich and 
oxygen-rich dust grains. There could also be an effect of ORD PNe originating from a range of pro- 
genitors, including very massive ones. One would expect to find PNe with low- to intermediate-mass 
progenitors in the upper right portion of the plot, where typically we find CRD PNe. 

In Figure 31 we plot the dust temperature against the nebular electron density. There is 
no particular segregation among the CRD and ORD or MCD PNe, apart from the high dust 
temperatures of CRD PNe as noted above. 

In Figure 32 we show the IR luminosity of the PNe vs. the nebular radius. If we could assume 
that the expansion velocity is uniform for all PNe, this plot would represent the evolution of IR 
luminosity in PNe. The expected decline of luminosity with radius is found as expected from a 
large PN population. Qualitatively, if the luminosity of the CS is mostly radiated through the 
dust continuum, and if the expansion is uniform and non accelerated, we can compare this plot 
with the log L/L Q vs. t [yr] plot in Stanghellini and Renzini (2000, Fig. 2 therein) and we see a 
correspondence. 

It is well known that nebulae and CS should be studied together to gain insight on their 
evolution and progenitor types. In the Galactic sample presented here the information about the 
CS is still very scarce. Given the compactness of the PNe in our sample, their CS magnitudes are 
typically unreachable from the ground, and the standard methods of temperature determinations 
from the Zanstra analysis are thus not applicable. Spectral types, on the other hand, are available 
for only very few CS of each dust class, thus correlations are not statistically sound. By correlating 
the CS spectral types (Weidmann & Gamen 2011) with the dust classes in this paper we noted that 
CS of the CRD PNe are of wels type in 6 of the 10 available spectra for this dust class, while most 
of the ORD PNe have Of or O(H) CS (8 of the 12 available spectra). While there is certainly a 
mild correlation between CS and dust class, the samples are too small to draw robust evolutionary 
scenarios. 

For a sizable part of the sample there will be high quality HST magnitudes available in the near 
future from the WFC3 images collected by us (Shaw et al., in preparation). For high excitation 
PNe, those whose CS are hot enough to doubly-ionize nebular helium, we can use the He II A4686 
flux as a probe of stellar temperature. In Table 4, column (5) we list the A4686 intensities found in 
the literature (see selection in Stanghellini & Haywood 2010), but the sample that includes reliable 
fluxes and dust parameters (temperature, IR luminosity) is too small for meaningful comparative 
analysis. 

We can use other criteria to characterize the hardness of the UV radiation of the CS both with 
IR and optical nebular emission lines. For example, Morgan (1984) showed that a good measure of 
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the hardness of the ionizing stellar flux in a PN was the excitation class, or EC, derived as EC=0.45 
1^5007/!/? for those PNe whose CS are not sufficiently hot to produce He II nebular emission, and 
EC=5.54 (Ia4686/I/9 +0.78) for the high excitation PNe. An estimate of the hardness of the UV 
field could be also obtained directly from the IRS spectra, for example by using the [Ne III] (15.56 
/im) to [Ne II] (12.81 /im) line fluxes ratio (Bernard-Salas et al. 2009). We list both the neon ratio 
described above, and calculated using our analysis of the IRS emission line spectra (Stanghellini 
et al., in prep.), and the excitation constant, calculated using the parameters in Stanghellini & 
Haywood (2010), in Table 4. Naturally neither of these criteria are ideal substitutes for a direct 
measurement of the stellar temperature, or even for an indirect measurement such as the stellar 
temperatures derived via the Zanstra method. For example, the neon flux ratio from IR emission 
lines is a reasonable stellar flux tracer only for intermediate excitation nebulae. Furthermore, the 
neon intensity ratio method fails if the PN is lumpy, non-homogeneous, or weather shocks dominate 
the emission. 

The [Ne III]/[Ne II] line ratio from the IRS spectra is a rough indication of the PN evolutionary 
stage, if the effects of shocks are absent or moderate. In fact, depending on where the line originate 
within the nebulae this ratio is tightly correlated to the excitation class or not, and the difference 
is expected to be larger in PNe with complicated shapes. In Figure 33 we show how the excitation 
class derived from the optical emission lines correlates quite well with the UV flux- sensitive neon 
intensity ratio. We observe in the Figure that the neon ratio is a better estimate of the PN 
excitation class in CRD and MCD PNe, where it correlates with the excitation constant with 
correlation coefficient (between the two logarithmic values) of 0.88 and 0.94 for the CRD and the 
MCD PNe respectively. On the other hand, the correlations is more scattered for ORD PNe, where 
the overall correlation coefficient between the two logarithmic values is 0.75, and at high excitation 
the correlation is very scattered. This finding complies with the hypothesis that in high excitation 
ORD PNe the neon ratio is not a good indicator of the hardness of the stellar flux. The explanation 
could lie in the asymmetric morphological type of a good fraction of ORD PNe, which was hinted 
in the Magellanic Cloud sample. 

The morphologies of our sample of Galactic PNe according to their dust type is illustrated in 
Table 5, where we list, for each dust class, the number of PNe in each of the major morphological 
classes, following the classification scheme by Stanghellini et al. (2002). Sequentially, we give in 
the Table the percentage of PNe in the separate categories of round (R), elliptical (E), bipolar core 
(BC), bipolar (B) and point-symmetric (P) PNe, as well as the cumulative classes of symmetric 
(round/elliptical) and asymmetric (bipolar /bipolar core/point-symmetric) PNe. It is clear that the 
few F PNe in this Table do not constitute a statistically significative sample. Among the other dust 
types there are some notable differences in the morphological distribution: CRD PNe are for the 
major part morphologically symmetric, with just a quarter of the sample displaying asymmetric 
features; on the other hand, both ORD and MCD PNe are almost equally distributed among the 
symmetric and asymmetric classes. It is also worth noting that the peculiar class of point-symmetric 
PNe are well populated by MCD targets, a very interesting coincidence of two peculiarities which 
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Fig. 31. — Tdust vs. Electron density, calculated from [S II] lines. Symbols are as in Figure 27. 
Temperature errorbars are as in Fig. 29, and for the electron densities it is assumed a 0.1 dex error 
(see text). 
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Fig. 32. — IR luminosity against nebular radius, in pc. Symbols are as in Figure 27. The error 
bar on the left lower corner includes the assumption of a typical distance uncertainty from the 
statistical distance derivation. 
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Fig. 33. — The PN excitation class, from optical emission lines (see text), versus the ratio of the 
IR line intensities of [Ne III] at A15.56 //m and [Ne II] at A12.81 /im . Symbols as in Figure 27. 
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is worth further study. 

5. Comparison of the Galactic and the Magellanic Cloud PN samples 

In the previous sections we have shown the atlas of the spectra of ^150 compact Galactic 
PNe that we plan to study in much more detail. In this data paper we aim to show the general 
characteristics of this sample, and compare it to the Magellanic Cloud PN IRS sample. Group 
and subgroup properties have been illustrated above; here we discuss what could be the leading 
evolutionary paths and the progenitor characteristics that might lead to the observed configuration. 
A thorough analysis of these issues is planned for a future paper. 

The first notable comparison between the Galactic and Magellanic Cloud samples is the sta- 
tistical distribution of dust types within the different environments. In Table 6 we give, for each 
main dust class, the PN frequency in the Galaxy (this paper, consisting in 150 PNe, 26 of which 
might be bulge PNe) and in the Magellanic Clouds, where both the S07 and the Bernard-Salas et 
al. (2009) samples have been included, for a total of 66 PNe in both Clouds. The major selection 
difference from the targets in the samples is the compactness of Galactic PNe, while the Magellanic 
Cloud PN have been selected on PN optical brightness. Given the tight correlation between physi- 
cal size and brightness, we assume the selection criteria will not affect the following discussion. Of 
the 150 Galactic PNe only 25 do not show notable molecular /dust features in their spectra, apart 
from a weak dust continuum, while in the Magellanic Cloud F PNe are ~41 percent of the sample. 
Sloan et al. 's (2008) showed that the number of evolved, dust-free stars increases with decreasing 
metallicity, broadly in agreement with our findings. 

Among the spectra with notable grain features, Galactic PNe are distributed almost uniformly 
among the CRD, ORD, and MCD dust classes. In the Magellanic Clouds, on the other hand, 
most PNe are CRD, only a few are ORD, and none have mixed-chemistry dust. In S07 we showed 
that the CRD PNe are also enriched in gas-phase carbon, and they represent the progeny of the 
intermediate mass in the PN progenitor range. The results of Table 6 indicate that PNe evolving 
from Magellanic Cloud progenitors are more likely to have carbon-rich dust than Galactic PNe, 
probably indicating that the Galactic disk PN progenitors are on average more massive than their 
homologous in the Clouds. Interestingly, none of the IRS observations of Magellanic Cloud PNe 
(including those by different Authors) show any indication of MCD class PNe, while this type is 
notable in the Galaxy. Whatever the mechanism to produce MCD PNe, it does not appear to be 
efficient at the Magellanic Cloud metallicity. 

The PNe in our Galactic sample that might belong to the bulge have characteristics that are 
even more extremely different from the Magellanic Clouds than those of Galactic disk PNe. In 
fact, nearly all bulge PNe have dust features, and only 11% of them have CRD types, while ORD 
(38%), plus especially MCD (46%) PNe are the majority. There seems to be continuity of dust 
class distribution from the bulge, to Galactic disk, to Magellanic Cloud PNe, where all the class 
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Fig. 34. — Dust temperature plotted against linear nebular dimensions (in cm). Triangles: F, 
diamonds: CRD. Filled symbols, like in Figures 27 through 35, represent Galactic PNe. Open 
symbols represent the Magellanic Cloud sample from S07. Solid lines: best fit to the CRD PN loci 
in the two samples. 
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Fig. 35. — Dust temperature vs. the integrated IR luminosity for Galactic and Magellanic PNe. 
Symbols are as in Figure 36. 
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Fig. 36. — Tdust vs. Electron density, calculated from [S II] lines, for Galactic and Magellanic PNe. 
Symbols are as in Figure 36. 
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frequencies change monotonically from one population to the next. It seems that metallicity plays 
a fundamental role in this picture: going from high (the bulge is a-element enriched with respect 
to the Galactic disk, Zoccali et al. 2006) to low (the Magellanic Clouds are half or quarter solar in 
a-element abundances) metallicity, the number of dust-rich PNe decreases, as does the fraction of 
PNe with ORD and MCD, while the number of CRD PNe increases. 

In Figure 34 we show the dust temperature vs. physical radii of the Galactic disk (filled 
symbols) and the Magellanic Cloud (open symbols) PNe, respectively from this paper (Galactic disk 
only) and S07. Given that we did not find MCD PNe in the Magellanic Clouds, and that the ORD 
PNe were a striking minority, here we only plot the F and CRD PNe for a meaningful comparison 
of the two samples. Distances, and consequently linear radii, are very sound for Magellanic Cloud 
PNe, being based on the independently determined distances to the Clouds, while in the Galactic 
case we use the statistical distance scale, which carries much larger, if reasonable, uncertainties. 
Nonetheless, the qualitative behavior or the CRD PNe in the two samples is similar in this plot, 
indicating a possible similarity of dust evolution in very different metallicity environments. By 
comparing the linear radii we see that we are sampling comparable expansion phases (assuming 
homogeneous expansion velocities), with the notable difference of few Galactic PNe of very small 
radii; this is sensible if we recall that the Galactic PN of the present sample were selected to be 
compact, unevolved ones. In spite of the uncertainties, and the caveats associated with adopting a 
simple, black-body temperature for the dust continuum, we note that in the correlations for CRD 
PNe are very good if considered separately for each host galaxy, with the temperatures in the 
Magellanic Cloud sample being higher than in the Galactic sample at the same radii. The slope of 
the correlation is slightly higher in the Magellanic Cloud (Tdust oc Rpn~ ' 32 ) than in the Galactic 
(Tdust of Rpn -0 ' 25 ) case, in both cases the correlation is tight, with correlation coefficient of -0.99 
and -0.75 for the Galactic disk and the Magellanic Cloud PNe respectively. If we assume that the 
expansion rates are uniform in two samples then we conclude that, at similar evolutionary stages, 
Magellanic PNe have higher dust temperatures than Galactic counterparts. If we assume that most 
CRD PNe have similar progenitor mass and metallicity (in each of the two samples shown here) 
we can conclude that the cooling of PN dust depends on the metallicity of the population. In 
particular, if the PN dust grains in the Magellanic Clouds are smaller than those in Galactic PNe 
due to metallicity, as seen in the ISM (e. g., Sandstrom et al. 2011), they are expected to retain 
higher dust temperatures (Li & Draine 2001), which is compatible with our analysis. 

It is worth mentioning that the models of stellar evolution (Marigo et al. 1999; Karakas et al. 
2003) predict that stars will experience the third dredge-up and the hot-bottom burning at lower 
masses for lower metallicities. There will therefore be a difference in the carbon-rich star progenitor 
masses for Galactic and Magellanic Cloud PNe. If the correlation between carbon abundances and 
carbon dust, demonstrated in Magellanic Cloud PNe by S07, will hold for the Galactic PNe as well, 
then the sequences of Galactic and Magellanic Cloud CRD PNe shown in Figure 34 could represent 
slightly different initial mass domains, which should be then taken into account in the comparison. 
To date, the sample of CRD Galactic PNe with a carbon abundance determination is way too small 
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for such an analysis to be meaningful. 

In Figure 35 we show the Td us t vs/ IR luminosity for the Galactic (filled symbols) and Magel- 
lanic Cloud (open symbols) PNe, for the F and CRD PNe. Both Galactic and extragalactic CRD 
PNe seem to populate a sequence where high luminosity correspond to high temperatures and vice 
versa. Galactic PNe have lower Tdust for similar luminosity, which we think is due to metallicity 
effect on dust grain behavior, as illustrated above, implying the existence of smaller carbonaceous 
grains at low metallicity. Figure 36 shows the relation of dust temperature with electron density, 
equally dense CRD PNe seem to have hotter dust in the Magellanic Cloud sample than the Galactic 
sample. 

6. Conclusions 

A large number of PNe with Spitzer/IRS spectra are presented in this paper. With the present 
sample, the early properties of ejected dust in Galactic planetary nebulae could be studied with 
sound statistical significance. The target sample includes all Galactic disk PNe with angular size 
smaller than ~4 arcsec. A few targets turned out not to be bona fide PNe, leaving us with a pure 
sample of 150 Galactic PNe to study. All PN spectra show dust continua and nebular emission 
lines. Solid state emission superimposed on the continuum displays a variety of characteristics. 
The dust classification, performed similarly to what we did for the Magellanic Cloud Spitzer/IRS 
spectra, consists of four major classes, determined by the type of dust that is prevalent in the 
spectra, including CRD, ORD, MCD, and featureless spectra PNe. Molecular /dust emission is 
much more apparent in this sample than in the Magellanic Cloud PNe, where featureless spectra 
with little or no continuum comprised a large fraction of the sample. 

The most populated dust class in this Galactic sample is the ORD (or oxygen-rich dust) class, 
including PNe featuring crystalline and amorphous silicates, and a few which show both types 
of grains. Second most populated class is the mixed-chemistry PNe, where both carbon-rich and 
oxygen-rich dust features are evident above the dust continuum. In contrast, in the Magellanic 
Clouds the ORD PNe were a tiny minority, and MCD PNe were absent. 

By separating the pure disk population and the possible bulge PNe within the present sample 
we see that the bulge PNe are mostly of dual chemistry, as previously observed. Bulge, disk, 
and Magellanic Cloud PNe form a sequence where decreasing metallicity determines a decreasing 
fraction of ORD and MCD PNe, and an increasing fraction of CRD PNe. 

Our analysis shows that CRD PNe define a rather narrow sequence when IR luminosity and 
physical radii are plotted against the dust temperature, which we interpret as an evolutionary 
sequence. The consequence of these findings is that CRD PNe must originate from progenitors 
with narrowly distributed mass, probably below ~3 M Q , or, non-type I PNe (Peimbert 1978). 
ORD and MCD PNe might not display a similar sequence because oxygen-rich grains radiate more 
efficiently than their CRD counterparts and their dust temperatures are lower than those of CRD 
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PNe. 

When comparing the Galactic and Magellanic Cloud CRD PN, in particular their T^ust evolu- 
tion and sequences, we find that dust temperatures are somewhat higher in the extragalactic than 
Galactic PNe for similar radii, which is in broad agreement with smaller carbonaceous dust grains 
observed at low metallicity in the ISM, and consequently possible smaller radiation efficiency in 
the Magellanic Cloud PNe than in their Galactic counterparts. Another likely explanation for the 
observed correlation is the effect of line blanketing at lower metallicity 

The data presented here have been already analyzed for the presence of some complex molecules, 
such as Ceo (fullerene, Garcia- Hernandez et al. 2010). More molecular analysis is performed at this 
time, as this represents the largest and most complete and homogeneous Spitzer/IRS Galactic PN 
data set available. We foresee that our group and others will take full advantage of the reduced 
data presented here for detailed future analysis of dust features. 
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i ni or p - /i 

19n35m54. 


,6s 


i o /i 1 r /( /i o 

+24d54m48s 


m n ,i nn 

010400 


1 


i o r o 

12.58 


TDAT 

FIN 


G 060.5+01.8 


TT O .1 A n 

He 2-440 


i ni o o no 

19n38m08. 


,5s 


i or J 1 c n^. 

+25dl5m40s 


nin/l nn 

010400 


-I 

1 


1 o r o 

12.58 


TDAT 

FIN 


/~*if\£iO o no o 

GU63.8-U3.3 


K 3-54 


oal n /i . a o 

2Un(J4m4o. 


,6s 


+ 25(101111,303 


n 1 n /i nn 

U1U4UU 


n 

u 


1 o r o 

12.58 


TDAT 

FIN 


GU67.9-UU.2 


T^ O CO 


2Un(Jomll. 


.OS 


i ') n .1 ') o ') i , 

+3Ud32m34s 


n 1 n /i nn 

U1U4UU 


1 


1 o r o 

12.58 


TDAT 

FIN 


G068.7+01.9 


T7" /I ^ 1 

K 4-41 


i ni r /"* o a 


,2s 


1 OOJOO 1 . 

+32d22ml2s 


nin/i nn 

010400 


l 


1 o r o 

12.58 


TDAT 
FIN 


GU68.7+14.8 


C~ A 1 

bp 4-1 


i ni^ nnw, nn 


,2s 


+38d2bmU7s 


m n 1 nn 

U1U4UU 


U 


44. U4 


TDAT 

FIN 


G069.2+02.8 


Ty' O Af\ 

K 3-49 


i ni r /i nn 

19n54m00 


,8s 


1 O O J OO 1 o 

+33d22ml2s 


niA/i nn 

010400 


1 


1 o r o 

12.58 


T~> AT 

PJN 


G077.7+03.1 


Kj 2 


om 1 r oo 

20nl5m22 


,4s 


i A n i o /( i /•> 

+40d34m46s 


nnn /( nn 

009400 


1 


12.58 


TDAT 

FIN 


G 079.9+00.4 


K 3-56 


oni c\n oi 

20n06m21 


."7- 

,7s 


i A A A AC\ 

+44dl4m40s 


nnn a nn 

009400 





1 O f O 

12.58 


DAT 

PJN 


G082.5+11.3 


NG 6833 


"i ni A c\ i . i 

19n49ml2. 


,5s 


+48d59m60s 


ai i a nn 

011400 





A A r\ A 

44.04 


T~i AT 

PJN 


G088.7+04.6 


TV O TO 

K 3-78 


oni /i r - 

20n44m45. 


,5s 


i f n 1 o o f o 

+50d23m52s 


nnn A nn 

009400 





12.58 


PN 


G095. 2+00.7 


K 3-62 


21h31m50 


,4s 


+52d33m52s 


009400 


1 


12.58 


PN 


G097.6-02.4 


M 2-50 


21h57m06. 


,9s 


+51d44m22s 


009400 





12.58 


PN 


G104.1+01.0 


Bl 2- 1 


22h20ml6. 


,7s 


+58dl4ml8s 


010600 


1 


12.58 


PN 


G107.4-00.6 


K 4-57 


22h48m34. 


,5s 


+58d29ml0s 


010600 


1 


12.58 


PN 


G107.4-02.6 


K 3-87 


22h54m23. 


,6s 


+56d43m38s 


009500 





12.58 


PN 


G112.5-00.1 


Kj 8 


23h23m22 


,8s 


+60d59ml8s 


009500 





44.04 


PN 


G184.0-02.1 


M 1- 5 


05h46m50. 


Is 


+24d22m01s 


010700 


1 


12.58 


PN 


G205.8-26.7 


MaC 2- 1 


05h03m39 


,7s 


-06d04m01s 


010700 





44.04 


PN 


G235.3-03.9 


M 1-12 


07hl9m21. 


,6s 


-21d43m57s 


010800 


1 


12.58 


PN 


G263.0-05.5 


PB 2 


08h20m52 


,6s 


-46dl7ml9s 


010900 





12.58 


PN 


G264.4-12.7 


He 2- 5 


07h47m21. 


,5s 


-51d09m01s 


010800 





12.58 


PN 


G274.1+02.5 


He 2- 34 


09h41ml3 


,9s 


-49d22m46s 


011400 


1 


12.58 


PN 


G275.3-04.7 


He 2- 21 


09hl4m04. 


,7s 


-55d22m28s 


011800 





44.04 


PN 


G278.6-06.7 


He 2- 26 


09hl9m44. 


,4s 


-59d06m22s 


011900 





12.58 


PN 


G285.4+01.5 


Pe 1- 1 


10h38m27 


,6s 


-56d47m05s 


009400 


1 


12.58 
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Table 1 — Continued 





Name 


alias 


G?9000 




09000 


IRS campaign 


mode 


t'cxp [s] 


(1) 






(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


PN 


G285. 


.4+02.2 


Pe 2- 7 


10h40m56. 


,6s 


-56dl4m24s 


009400 





12. 


.58 


PN 


G286 


.0-06.5 


He 2- 41 


10h07m58. 


9s 


-63d49m52s 


011000 





12. 


.58 


PN 


G289 


.8+07.7 


He 2- 63 


Ilh24ml7. 


Is 


-52d45m47s 


011000 





44. 


.04 


PN 


G294. 


.9-04.3 


He 2- 68 


Ilh31m45. 


9s 


-65d58ml5s 


011100 


3 


12. 


.58 


PN 


G295 


.3-09.3 


He 2- 62 


Ilhl7ml4. 


,2s 


-70d55m06s 


009400 





12. 


.58 


PN 


G296 


,3-03.0 


He 2- 73 


Ilh48m38. 


,7s 


-65d08m39s 


011100 


3 


12. 


.58 


PN 


G297 


.4+03.7 


He 2- 78 


12h09m24, 


,5s 


-58d48m23s 


011400 





12. 


.58 


PN 


G300 


.7-02.0 


He 2- 86 


12h30m30. 


,5s 


-64d52m07s 


012000 


1 


12. 


.58 


PN 


G307 


.5-04.9 


MyCn 18 


13h39m34, 


9s 


-67d22m50s 


009400 


1 


12. 


.58 


PN 


G309 


.0+00.8 


He 2- 96 


13h42m36. 


0s 


-61d22m28s 


009400 


1 


12. 


.58 


PN 


G309 


.5-02.9 


MaC 1- 2 


13h54m37. 


,5s 


-65d05m33s 


009400 





44. 


.04 


PN 


G311 


.1+03.4 


He 2-101 


13h54m55. 


,6s 


-58d27ml5s 


009400 


1 


12. 


.58 


PN 


G321 


.3+02.8 


He 2-115 


15h05ml6. 


,7s 


-55dllm09s 


010600 


1 


12. 


.58 


PN 


G324 


.2+02.5 


He 2-125 


15h23m33. 


,8s 


-53d57m30s 


010600 





12. 


.58 


PN 


G324 


.8-01.1 


He 2-133 


15h41m58. 


,7s 


-56d36m24s 


010600 


1 


12. 


.58 


PN 


G325 


.0+03.2 


He 2-129 


15h25m41. 


,8s 


-52d56m32s 


009500 





12. 


.58 


PN 


G325 


.8-12.8 


He 2-182 


16h54m35. 


0s 


-64dl4m27s 


010600 


1 


12. 


.58 


PN 


G326 


.0-06.5 


He 2-151 


16hl5m42. 


Is 


-59d53m60s 


010600 


1 


12. 


.58 


PN 


G327. 


.1-01.8 


He 2-140 


15h58m08. 


0s 


-55d41m49s 


010600 


1 


12. 


.58 


PN 


G327 


.8-06.1 


He 2-158 


16h23m36. 


Is 


-58d25m23s 


010600 





12. 


.58 


PN 


G327. 


.9-04.3 


He 2-147 


16hl4m00. 


9s 


-56d59m26s 


010600 


1 


12. 


.58 


PN 


G329 


.4-02.7 


He 2-149 


16hl4m27. 


,7s 


-54d53m41s 


010600 





12. 


.58 


PN 


G331 


.0-02.7 


He 2-157 


16h22ml4. 


Is 


-53d40m53s 


010600 


1 


12. 


.58 


PN 


G334 


.8-07.4 


SaSt 2- 12 


17h03ml3. 


,2s 


-54d01m46s 


010600 





12. 


.58 


PN 


G336 


.3-05.6 


He 2-186 


16h59m45. 


,2s 


-51d47m59s 


010600 





12. 


.58 


PN 


G336 


.9+08.3 


StWr 4-10 


16h02ml4. 


,5s 


-41d39m39s 


010600 





12. 


.58 


PN 


G340 


.9-04.6 


Sa 1-5 


17hllm36. 


,7s 


-47d30m52s 


010600 





12. 


.58 


PN 


G341 


.5-09.1 


He 2-248 


17h36ml6. 


9s 


-49d31m35s 


010600 





44. 


.04 


PN 


G343 


.4+11.9 


H 1- 1 


16hl3m30. 


6s 


-34d41m41s 


010600 





44. 


.04 


PN 


G344 


.2+04.7 


Vd 1-1 


16h42m38. 


,4s 


-39d00m30s 


010600 





12. 


.58 


PN 


G344 


.4-06.1 


Wray 16-278 


17h30ml5. 


0s 


-45d28m34s 


010600 





44. 


.04 
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Table 1 — Continued 



Name 


alias 


«2000 


<52000 


IRS campaign 


mode 


t'cxp [s] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 





O A A 


A 1 no Q 

.4+02.8 


A T A 1 E 

Vd 1-5 


lbn51mo7. 


o . 
OS 


-4UaU5mobs 


ai n/?nn 

UlUbUU 


n 
(J 


A A 

44 


n a 
.1)4 


T~i AT 

PJN 


G344 


,8+03.4 


Vd 1-3 


1 c 1 /in A~\ 

16n49m41. 


9s 


on 1 o / ' p" 

-39d26m57s 


nnn p* nn 

009500 





A i 

44. 


n a 

.04 


T~i AT 

PIN 


G345 


.0+1)4.3 


^ J A 1 O 

Vd 1-2 


1 Ac r n 

16n46m50. 


OS 


O O J /( O rr r" 

-38d42m55s 


ni n^nn 

010600 





1 O 

12. 


CT O 

.58 


T~i AT 

PIN 


G348 


A r\ A i 

.4-04.1 


H 1-21 


1 ^71 O O CT T 

17n32m57. 


,7s 


A A J\f\ A 1 /< ._. 

-41d04ml4s 


ni c\rr\r\ 

010600 





1 O 

12. 


CT O 

.58 


T~> AT 

PJN 


G348 


o nn n 

.8-09.0 


TT o on/"* 

He 2-306 


17n56m46. 


2s 


/i o l no rr 

-43d08m55s 


ni npan 

010600 





1 O 

12. 


.58 


T~i AT 

PIN 


G350 


o no a 

.8-02.4 


TT 1 on 

H 1-22 


1 ^71 OO OO 

17n32m22. 


n 

0s 


OTJrf7 00„ 

-37d57m23s 


010600 


1 


1 O 

12. 


.58 


T~i AT 

FIN 


G351 


.3+07.6 


TT 1 A 

H 1- 4 


1 Pi CT O /II 

16n53m41. 


,4s 


01 J /ip OO 

-31d46m32s 


ni n^nn 

010600 





A A 

44. 


n a 

.04 




G351 


.y-oi. y 


Wray Lb-zob 


1 tyv, o o n n 

17noomUU. 


OS 


O Ct A A O , ^ r 1 „ 


ni n/^nn 

010600 


1 


1 O 

12. 


.58 




G352 


a i no n 

.6+03.0 


I T -1 O 

nl-o 


17nl4m45, 


,7s 


o o A on i o.. 


ni ncnn 

010600 


n 




A A 

44. 


.04 


TDAT 


G354. 


.2+04.3 


a /r o in 


17nl4ml<l 


,4s 


-31d25m3os 


ni ncnn 

010600 


n 




1 o 

12. 


r o 

.58 


TDAT 

PIN 


G354 


.9+03.5 


mi q ^ 

In 3- 6 


1 ti in on 

17nl9m20. 


,1s 


o i J i o on 

-31dl2m39s 


OlObOO 


l 


1 o 

12. 


r o 

.58 


TDAT 

PIN 


G355 


o no c 

,2-02.5 


tt i on 

H 1-29 


1 'TT A A 1 O 

17n44ml3. 


,7s 


-34dl7m32s 


ni n*° nn 

010600 


l 


1 O 

12. 


.58 


TDAT 

PIN 


G355 


.7-03.0 


TT 1 O O 

H 1-33 


1 Tl AT CT O 

17n47m58. 


,4s 


O /( J 1 O CT 1 

-34dl3m51s 


n i n/"*nn 

010600 





A A 

44. 


n a 

.04 


TDAT 

PIN 


Good 


.9+02.7 


in 3-10 


1 ti o /i /in 

17n24m40. 


,8s 


on J CT 1 CT O-. 

-30d51m58s 


n 1 npnn 

010600 


l 


1 O 

12. 


CT O 

.58 


T~> AT 

PJN 


G356 


.2-04.4 


Gn 2-1 


1 ^71 C ^ oo 

17n54m32, 


9s 


o A 1 oo on 

-34d22m20s 


n i r\r* nn 

010600 


l 


12. 


.58 


TDAT 

PIN 


Good 


.5-03.0 


TT O 0^7 

H 2-27 


1 ^Tl CT 1 CT n 

17n51m50. 


,5s 


OO J yl7 OT„ 

-33d47m35s 


n i n*° nn 

010600 


l 


1 O 

12. 


.58 


PN 


G356 


.5+01.5 


Th 3-55 


17h30m58. 


,7s 


-31d01m05s 


010600 


l 


12. 


.58 


PN 


G356 


.8+03.3 


Th 3-12 


17h25m06, 


,0s 


-29d45ml6s 


010600 


l 


12. 


.58 


PN 


G357. 


.1-06.1 


M 3-50 


18h04ml5. 


,4s 


-34d34ml8s 


010600 





44. 


.04 


PN 


G357 


.1+01.9 


Th 3-24 


17h30m57. 


,9s 


-30d23m05s 


010600 





44. 


.04 


PN 


G357. 


.2+02.0 


H 2-13 


17h31ml4. 


,6s 


-30dl6m21s 


010600 





44. 


.04 


PN 


G357 


.6+01.7 


H 1-23 


17h32m47. 


,0s 


-30d00ml6s 


010600 


1 


12. 


.58 


PN 


G357 


.6+02.6 


H 1-18 


17h29m42. 


,7s 


-29d32m49s 


010600 


1 


12. 


.58 


PN 


G358 


.2+03.6 


M 3-10 


17h27m20. 


,1s 


-28d27m50s 


010600 


1 


12. 


.58 


PN 


G358 


.3+01.2 


Bl B 


17h37m06. 


,7s 


-29d46m01s 


010600 





44. 


.04 


PN 


G358 


.5-04.2 


H 1-46 


17h59m02, 


,4s 


-32d21m42s 


010600 


1 


12. 


.58 


PN 


G358 


.5+02.9 


Al 2-F 


17h30m36, 


,8s 


-28d41m48s 


010600 





44. 


.04 


PN 


G358 


.6+01.8 


M 4- 6 


17h35ml3. 


,9s 


-29d03m09s 


010600 


1 


12. 


.58 


PN 


G358. 


.7-02.7 


Al 2-R 


17h53m45. 


,1s 


-31d31mlls 


010600 





44. 


.04 


PN 


G358 


.7+05.2 


M 3-40 


17h22m28. 


,2s 


-27d08m41s 


010600 


1 


12. 


.58 


PN 


G358. 


.9-03.7 


H 1-44 


17h58ml0. 


,5s 


-31d42m55s 


010600 


1 


12. 


.58 
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Table 1 — Continued 



Name alias 02000 ^2000 IPS campaign mode t exp [s] 

(1) (2) (3) (4) (5) ' (6) (7) 



PN G358.9+03.4 H 1-19 17h30m02.5s -27d59ml6s 010600 1 12.58 

PN G359.3+03.6 Al 2-E 17h30ml4.3s -27d30ml8s 010600 1 12.58 

PN G359.4+02.3 Th 3-32 17h35m22.1s -28dl2m54s 010600 44.04 
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Table 2. Classification scheme 



class 


description 


subclass 


N PN 


(1) 


(2) 


(3) 


(4) 


F 


dust continuum 




25 (17%) 


CRD 


carbon-rich dust 




38 (25%) 






aromatic 


13 






aliphatic 


22 






aromatic / aliphatic 


3 


ORD 


oxygen-rich dust 




45 (30%) 






crystalline 


16 






amorphous 


24 






crystalline / amorphous 


1 



MCD mixed-chemistry dust 



42 (28%) 



Table 3. Dust classification and dust parameters 



Name class sub. fit a F re f. Tjust log Ltr IRE 

a c [Jy] [K] [l ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



PN G000.3-02.8 


F 





B 


1.02 




65.81 


± 


2.58 






PN G000.6-01.3 


CRD 


3 


N 
















PN G000.6-02.3 


MCD 


7 


N 
















PN G000.8-01.5 


MCD 


7 


A 


2.20 


4.77 e 


79.38 


± 


1.80 






PN G000. 8-07.6 


CRD 


-1 h 


N 
















PN G001. 7-01.6 


MCD 


7 


B 


2.58 




76.23 


± 


5.56 






PN G001. 7-04.6 


ORD 


4 


A 


2.22 


1.41 d 


82.31 


± 


1.76 


2.920 


2.920 


PN G002.4-03.7 


MCD 


7 


A 


2.79 


5.21 e 


77.05 


± 


0.90 


2.837 


2.837 


PN G002.6-03.4 


MCD 


7 


A 


2.68 


9.19 e 


75.74 


± 


0.88 






PN G002.9-03.9 


ORD 


5" 


N 
















PN G003.0-02.6 


CRD 


2" 


N 
















PN G003. 1+03.4 


MCD 


7 


A 


1.95 e 


3.92 ± 0.49 d 


79.12 


± 


1.83 


3.064 


3.064 


PN G003. 2-04.4 


ORD 


5 b 


B 


2.39 




77.88 




3.31 


2.346 


2.346 


PN G003.9-14.9 


ORD 


r-h 


A 


5.52 


1.93 


70.27 




2.00 


2.714 


2.714 


PN G004.1-03.8 


ORD 


5 b 


B 


0.00 




126.00 


± 


4.55 


1.768 


1.768 


PN G004.3-02.6 


MCD 


7 


N 
















PN G004.9-04.9 


MCD 


7 


C 


0.16 


5.19 


98.87 


± 


3.18 


2.669 


2.669 


PN G005.5+0.27 


MCD 


7 


C 


6.63 


13.30 6 


57.07 


± 


1.13 






PN G006.0+02.8 


MCD 


7 


A 


6.38 


l.ll d 


58.13 


± 


0.74 






PN G006.3+04.4 


ORD 


5 b 


N 
















PN G006.8+02.0 


MCD 


7 


N 
















PN G006.8-19.8 


CRD 


3 


N 
















PN G007.5+04.3 


ORD 


4 


N 
















PN G008. 1-04.7 


ORD 


4 


N 
















PN G008. 2-04.8 


F 





N 
















PN G008.6-02.6 


ORD 


5 


N 
















PN G008.6-07.0 


F 





N 
















PN G009.3+04.1 


ORD 


6 b 


C 


2.19 


0.85 e 


87.89 


± 


1.40 






PN G010.6+03.2 


CRD 


1 


A 


1.58 


l.82 d 


91.72 


± 


2.19 






PN G011. 1+07.0 


ORD 


5 


N 

















Table 3 — Continued 



Name class sub. fit a F re f. Td us t log Lir IRE 

a c [Jy] [K] [l ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



PN 


G011. 3-09.4 


ORD 


•5 


N 






PN 


G011.7-06.6 


NA 


8 


N 






PN 


G012.5-09.8 


CRD 


2 


N 






PN 


G012.6-02.7 


MCD 


7 


N 






T~i AT 

PJN 


G014.0-05.5 


1 





C 


1 


34 


DAT 

PJN 


G014.3-05.5 


CRD 


1 


AT 

JN 






T~i AT 

PJN 


G016. 9-02.0 


TV I / " ^ I ~\ 

MCD 


7 


AT 

JN 






T~i AT 

PJN 


G018.6-02.2 


1 





AT 

JN 






DAT 

PJN 


G019.2-02.2 


/ \ T ) TA 

ORD 


r* 






3 


28 


T~i AT 

PJN 


G019.4-05.3 


TV T/^lTA 

MCD 


7 


AT 

JN 






PN 


G019.7+03.2 


TV T/"lTA 

MCD 


7 


A 

A 


1 


00 


T~» AT 

PJN 


G019.8+05.6 


1 





AT 

JN 






PN 


G023.8-01.7 


MCD 


7 


A 

A 


3 


36 


T~» AT 

PJN 


G023.9+01.2 


TV T T A 

MCD 


7 


A 


2 


85 


PN 


G025. 3-04.6 


ORD 


5 


N 






PN 


G026.0-01.8 


ORD 


5 fe 


N 






PN 


G027.6-09.6 


ORD 


5 


N 






PN 


G031.0+04.1 


ORD 


4 


A 


4 


12 


PN 


G032.5-03.2 


MCD 


7 


A 





88 


PN 


G032.9-02.8 


CRD 


2 


A 





00 


PN 


G034.0+02.2 


ORD 


4 


A 





92 


PN 


G038.4-03.3 


F 





A 


1 


62 


PN 


G038.7-03.3 


F 





A 


2 


11 


PN 


G041. 8+04.4 


CRD 


3 


A 


1 


50 


PN 


G042.0+05.4 


ORD 


5 


N 






PN 


G042.9-06.9 


ORD 


6 


A 


7 


16 


PN 


G044.1+05.8 


NA 


8 


N 






PN 


G045.9-01.9 


MCD 


7 


A 


2 


98 


PN 


G048.1+01.1 


CRD 


1 


N 






PN 


G048.5+04.2 


CRD 


2 fe 


B 





00 



1.91 


± 


0.58 d 


77.46 


± 


3.42 












1 0.Oo 


± 


A t^A 
4.04 


z.yio 


Q1 K 

z.yio 


LU.o ( 


it 


z.lo 


i no /in 
1UZ.4U 


i 

± 


O.OD 


Q O/IO 
O.Z4Z 


Q OAO 


8.77 


zb 


2.11 d 


74.76 


± 


1.83 


3.377 


3.377 


10.93 


± 


1.70 d 


76.23 


± 


1.49 






5.50 


± 


0.54 d 


85.20 


± 


2.53 


3.218 


3.218 






1.36 d 


104.20 


± 


6.46 










2.56 e 


146.80 


± 


7.62 


3.000 


3.000 






3.56 e 


88.48 


± 


1.87 


2.699 


2.699 


2.14 


± 


0.01 d 


84.15 


± 


3.54 






3.23 


± 


0.73 d 


80.52 


± 


3.25 










0.52 d 


127.20 


± 


5.24 










0.65 d 


63.87 


± 


0.66 


3.703 


3.703 






3.46 e 


77.65 


± 


0.79 


3.394 


3.394 








108.20 


± 


3.60 


2.153 


2.153 



Table 3 — Continued 



Name class sub. fit a F re f. Tjust log Lir IRE 

a c [Jy] [K] [l ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



pw pn^l D-l_n9 8 
.TIN V-rUOl.U-rUZ.o 




A 

4 


A 


n nn 


o.oy it u.04: 


994 Rn -i- 1 9 e;q, 
Zi4. ou it iz.oy 






pw pn c i9 Q_i_n9 7 

.rlN vjUOZ.y^UZ. I 


PR Pi 


O 

Z 


M 
in 












pm pnw q no 7 

.r IN vjUOZ.y-UZ. / 


P 


n 
u 


M 












pm pn^ 1 ni 8 

rlN VjUOO.l-Ul.o 


p 


n 
u 


M 
IN 












.TIN vjUOO.O-UU.O 


PRPi 


o 
z 


M 
IN 












pm Pn^s Q^m ^ 
-TIN vjUOo.y^Ul.o 


DP Pi 


A 

4 


M 

IN 












.rlN vjUOy .4:^UZ.o 


PRD 


9 
Z 


A 
f\ 


n nn 
u.uu 


n 8n<^ 

U.oU 


1 9/i nn 4- A 8Q 
iz4.uu it 4:.oy 


9^19 
Z.01Z 


9 M 9 
Z.OIZ 


pm pn^Q Q-i-n9 n 

-TIN VjUOy.y^UZ.U 


PRD 

V III/ 


9 

z 


^N 












-TIN vjUOU.O^UI.O 


PR Pi 


/I 


A 
. \ 


9 ^1 
Z.ol 


9 7n -i- n 1 

Z. / U It U. 1U 


oc Q/i _i_ 9 ni 
otj.y^J: it z.ui 


9 9Q7 

o.zy ( 


1 9Q7 

o.zy / 


.TIN vjUUO.O-UO.O 


PRD 
v 1 1 1 / 


9 
z 


IN 












pm pdr7 q nn 9 
.tin v^uoi.y-uu.z 


PRD 




P 


n nn 
u.uu 


oo.y4 it z. i u 


1 0A KC\ -4- 1 03 

iz4.ou it i.yo 




0. 104 


pm Pnfis 7^m o 

-TIN viUOo.f+Ul.y 


PR Pi 


,1 


A 
A 


1 .40 


Z.4:4: It U.1U 


o« o« _|_ i 71 
oO.ZO it 1. / 1 


9 an^; 

Z.oUO 


9 sn^ 

Z.oUO 


PM PnfiR 7-Ll/l 8 


PRPi 


1 


M 
IN 












PM fflfiQ 9^n9 8 
.TIN V^rUDy.Z-j-UZ.O 


PR Pi 





A 
A 


^ n^ 


l.±U it U.oO 


Ol 98 _|_ 9 on 

oi.zo it z.oy 


0.001 


0.001 


"oat r^r\ r i'i 'i i no 1 

FIN G077. 7+03.1 


1 / 





IN 












PN G079.9+06.4 


ORD 


5 b 


B 


0.70 




118.00 ± 4.39 






PN G082.5+11.3 


F 





N 












PN G088.7+04.6 


F 





A 


1.59 


1.99 ± 0.36 d 


84.24 ± 3.63 


2.550 


2.550 


PN G095.2+00.7 


CRD 


2 


N 












PN G097.6-02.4 


ORD 


5 b 


N 












PN G104.1+01.0 


CRD 


1 


A 


0.00 


2.66 ± 0.10 d 


157.70 ± 3.83 


3.114 


3.114 


PN G107.4-00.6 


NA 


8 


N 












PN G107.4-02.6 


CRD 


t 


N 












PN G112. 5-00.1 


F 





N 












PN G184.0-02.1 


CRD 


2 


N 












PN G205.8-26.7 


CRD 


2 


B 


1.02 




122.30 ± 13.89 






PN G235.3-03.9 


CRD 


2 


N 












PN G263.0-05.5 


CRD 


2 


N 












PN G264.4-12.7 


CRD 


2 


N 












PN G274.1+02.5 


NA 


8 


N 
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Table 3 — Continued 



Name class sub. fit a F re f. Tdust log Lir IRE 

a c [Jy] [K] [l ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



T~i AT 

PJN 


G275 


,3-04.7 


CRD 


2 


AT 

JN 






T~i AT 

PJN 


/""I OTO 

G278. 


,6-06.7 


/"1T> TA 

CRD 


2 


A 

A 


0, 


.53 


PN 


G285. 


.4+01.5 


/~1T"» TA 

CRD 


1 


N 






T"» AT 

PJN 


G285. 


,4+02.2 


1 





AT 

JN 






T~» AT 

PJN 


r \ nor' 

G286. 


.0-06.5 


^~1T» TA 

CRD 


2 


A 

A 


0, 


,00 


T~i AT 

PJN 


G289 


,8+07.7 


1 





13 


0, 


,00 


T~i AT 

PJN 


G294 


.9-04.3 


^~1T» TA 

CRD 


2 


AT 

JN 






T~i AT 

PJN 


G295 


.3-09.3 


AT) TA 

ORD 


5 


A 

A 


1. 


TO 

.78 


T~i AT 

PJN 


G296 


.3-03.0 


MCD 


7 


A 

A 


2 


OT 

.27 


T~i AT 

PJN 


G297 


A i C\ O T 

.4+03.7 


^1T» TA 

CRD 


1 


AT 

JN 






PN 


G300 


,7-02.0 


TV I I ~\ 

MCD 


7 


A 

A 


3 


.65 


PN 


G307. 


.5-04.9 


MCD 


7 


A 

A 


2 


.42 


PN 


G309 


.0+00.8 


MCD 


7 


C 


2 


.07 


T~i AT 

PJN 


G309 


.5-02.9 


^~1T> TA 

CRD 


i b 
1 


AT 

JN 






PN 


G311 


.1+03.4 


NA 


8 


N 






PN 


G321 


.3+02.8 


CRD 


2 


N 






PN 


G324 


.2+02.5 


MCD 


7 


A 


3 


.26 


PN 


G324 


.8-01.1 


MCD 


7 


A 


3 


.12 


PN 


G325 


.0+03.2 


ORD 


5 


A 





.98 


PN 


G325 


.8-12.8 


ORD 


5 


A 


5 


.63 


PN 


G326 


,0-06.5 


ORD 


5 


A 


7. 


,46 


PN 


G327 


,1-01.8 


MCD 


7 


A 


3 


.12 


PN 


G327 


,8-06.1 


MCD 


7 


A 


0. 


,68 


PN 


G327 


,9-04.3 


NA 


8 


N 






PN 


G329 


.4-02.7 


F 





N 






PN 


G331 


.0-02.7 


ORD 


4 


A 


3 


,26 


PN 


G334 


.8-07.4 


ORD 


4 


A 


0. 


,81 


PN 


G336 


.3-05.6 


CRD 


1 


N 






PN 


G336 


.9+08.3 


ORD 


5 


N 






PN 


G340 


,9-04.6 


F 





N 







1.52 


128.10 ± 13.15 


2.815 


2.815 


r\ f TP 

0.57 


143.80 


± 


8.65 


o Tr r\ 

2.750 


o Tr r\ 

2.750 




liO.OU 






1 T A O 

1.743 


1 T A O 

1.743 


r\ o Tf 5 

0.87 


97.33 


1 

± 


7.08 






5.00 ± 0.43 d 


82.11 


± 


2.26 


3.063 


3.063 


18.48 ± l.ll d 


71.93 


± 


1.62 


3.450 


3.450 


17.17 ± 0.78 d 


87.02 


± 


1.20 


2.834 


2.834 


18.27 ± 1.43 d 


90.02 


± 


2.37 






2.86 ± 0.56 d 


70.53 


± 


4.14 


3.167 


3.167 


19.45 ± 2.31 


76.59 


± 


1.67 


3.125 


3.125 


1.33 d 


99.32 


± 


6.09 


2.868 


2.868 


0.72 d 


65.55 


± 


0.59 


2.496 


2.496 


4.70 ± 0.93 d 


57.45 


± 


0.39 


3.560 


3.560 


11.77 ± 3.41 d 


73.83 


± 


1.28 


3.318 


3.318 


1.12 d 


103.90 


± 


6.45 


3.022 


3.022 


2.61 d 


75.65 


± 


1.59 


3.068 


3.068 



0.41 e 133.10 ± 10.75 



Table 3 — Continued 



Name class sub. fit a F re f. Tdust log Lir IRE 

a c [Jy] [K] [l ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



PN G341.5-09.1 


F 





N 
















TiAT /"I O A O /( i 11 r» 

PJN G343.4+11.9 


ORD 


5 


B 


1.28 




97.56 ± 14.71 


2.301 


2.301 


PW G344.2+04.7 


ORD 


5 


A T 

JN 
















PN G344.4+02.8 


CRD 


2 


AT 

N 
















PM G344.8+03.4 


/^T"» TA 

CRD 


1 

1 


AT 

JN 
















PIN G345.0+04.3 


URD 


5 


A 

A 


4.62 


n t A(l 

0.74 




4- 


"\ 9Q 






PIN G348. 4-04.1 


TV T f \ T "\ 

MOD 


7 


AT 

JN 
















t ~» \ t f \ o a o o r\rv n 

PIN G348. 8-09.0 


t 





A 

A 


1 ro 

1.53 


2.16 


OO.UO 


4 

ZIZ 


9 1 7 
i. J i 






PIN G350.8-02.4 


URD 


4 


A 

A 


2.74 


1.60 


79 97 


4 

ZIZ 


1 9A 






PIN G351.3+07.6 


URD 


5 


AT 

IN 
















PN G351.9-01.9 


TV T f \ T A 

MCD 


7 


A" 

JN 
















PN G352. 6+03.0 


IV T f~ \ T A 

MCD 


7 


B 


2.38 




77 nn 
/ / .uu 


4 


^ 9^ 

o.zo 


2.988 


2.988 


PN G354.2+04.3 


CRD 


1 


AT 

JN 
















PN G354.9+03.5 


IV T T A 

MCD 


7 


AT 

JN 
















PN G355.2-02.5 


MCD 


7 


A 


0.22 


1.53 e 


124.30 


± 


4.33 


2.059 


2.059 


PN G355.7-03.0 


CRD 


1 


B 


1.08 




97.23 


± 


6.07 


2.770 


2.770 


PN G355.9+02.7 


MCD 


7 


C 


0.70 


7.81 ± 2.33 d 


96.69 


± 


5.66 


3.179 


3.179 


PN G356. 2-04.4 


CRD 


6 


A 


4.17 


2.98 ± 0.14 d 


74.39 




1.81 


3.062 


3.062 


PN G356.5+01.5 


ORD 


4 


C 


1.28 


7.21 e 


88.89 


± 


1.78 


2.446 


2.446 


PN G356.5-03.6 


F 





C 


0.46 


2.03 rf 


92.95 


± 


1.27 






PN G356.8+03.3 


MCD 


7 


N 
















PN G357.1+01.9 


F 





N 
















PN G357.1-06.1 


F 





N 
















PN G357.2+02.0 


MCD 


7 b 


B 


0.19 




118.70 


± 


7.02 






PN G357.6+01.7 


ORD 


4 


A 


2.29 


5.35 d 


83.75 


± 


2.01 


3.153 


3.153 


PN G357.6+02.6 


MCD 


7 


A 


1.86 


6.59 e 


86.34 


± 


2.49 


2.026 


2.026 


PN G358.2+03.6 


ORD 


4 


A 


2.08 


2.65 d 


89.06 


± 


2.60 


3.116 


3.116 


PN G358.3+01.2 


MCD 


7 


N 
















PN G358.5+02.9 


F 





B 


0.97 




110.80 


± 


5.39 






PN G358.5-04.2 


ORD 


6 


A 


3.86 


2.17 ± 0.35 d 


75.65 


± 


1.61 


1.616 


1.616 



Table 3 — Continued 



Name class sub. fit a F re f. T^ust log Lir IRE 

a c [Jy] [K] [L ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



PN G358.6+01.8 


MCD 


7 


A 


1.91 e 


4.28 


90.02 ± 2.42 


PN G358.7+05.2 


MCD 


7 


C 






73.37 ± 1.94 


PN G358.7-02.7 


NA 


8 


N 








PN G358.9+03.4 


MCD 


7 


A 


2.52 


5.17 ± 1.02 d 


79.98 ± 1.93 


PN G358.9-03.7 


MCD 


7 


N 








PN G359.3+03.6 


F 





A 


0.99 


1.23 d 


86.01 ± 1.79 


PN G359.4+02.3 


MCD 


7 


A 


2.14 


3.98 d 


85.92 ± 3.75 



a Subclasses: 0=featureless, l=aromatic, 2=aliphatic, 3=aromatic/aliphatic, 4=crystalline, 
5=amorphous, 6= crystalline/amorphous, 7=MCD, 8=other; b : Class or subclass uncertain. c 
Black body fit type: NC= not converging (see text); NA: not available (not a PN); A=good fit, 
excellent F re f (60 or 65 /an) constraints; B=good fit, no constraints available; C=good fit, poor 
constraints (not used in the plots). d F re f=F65^ m from Akari. e F m f=Fso lim from IRAS. * Flux 
from out grey-body fit, calculated either at 60 or 65 /im, depending on F re f. 
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Table 4. Other PN parameters 











J-4686 


TsJp vfifin 


EC 




fknrl 


PTT1 
|C111J 


[1,111 j 


[1/9—1001 






(1) 


(2) 




(41 


(5) 


(6) 


(7) 


PN G000.3-02.8 






2.511 




1.257 


4.181 


PN G000.6-01.3 






3.357 




-1.204 




PN G000.6-02.3 






3.280 








PN G000.8-01.5 






4.166 








PN G000.8-07.6 






3.304 




1.162 


4.415 


PN G001.7-01.6 






3.767 








PN G001. 7-04.6 


4.731 


17.46 


3.787 


2.5 


1.133 


2.839 


PN G002.4-03.7 


1.439 


17.39 


4.352 








PN G002.6-03.4 






3.722 








PN G002.9-03.9 








35.0 




6.260 


PN G003.0-02.6 


10.880 


17.63 




34.0 




6.205 


PN G003.1+03.4 


2.868 


17.47 


3.215 








PN G003. 2-04.4 


8.955 


17.55 






1.516 


4.127 


PN G003. 9-14.9 


0.547 


17.39 


3.722 




1.801 


3.740 


PN G004.1-03.8 


7.424 


17.57 






1.097 


4.095 


PN G004.3-02.6 


5.251 


17.71 






0.460 


2.209 


PN G004.9-04.9 


2.521 


17.50 


3.428 




-0.548 


0.112 


PN G005.5+0.27 












0.814 


PN G006.0+02.8 






2.707 








PN G006. 3+04.4 


2.515 


17.47 




4.5 


1.829 


6.422 


PN G006.8+02.0 










0.292 


1.696 


PN G006.8-19.8 






3.206 


13.4 




5.064 


PN G007.5+04.3 












1.255 


PN G008. 1-04.7 


4.401 


17.47 


3.240 




0.733 


2.416 


PN G008. 2-04.8 


1.967 


17.46 


3.240 




1.516 


4.181 


PN G008. 6-02.6 










1.614 


4.477 


PN G008.6-07.0 


7.938 


17.55 


2.825 


15.0 


1.264 


5.152 


PN G009.3+04.1 










1.535 


4.851 


PN G010.6+03.2 






4.114 




0.316 


1.561 


PN G011. 1+07.0 








50.0 




7.091 



Table 4 — Continued 



Name Rg log Rpn log N e I4686 Ne ratio EC 

[kpc] [cm] [cm" 3 ] [1/3=100] 

(1) (2) (3) " (4) (5) (6) (7) 



PN G011.3-09.4 


1.822 


16.98 


4.114 




-1.220 


0.283 


PN G011.7-06.6 


4.816 


17.79 


4.270 








PN G012.5-09.8 


2.959 


17.45 




7.8 




5.845 


PM G012.6-02.7 


3.112 


17.59 


A nro 

4.058 








PM G014.0-05.5 








on 

29.0 


1.671 


5.928 


PIN G014.3-05.5 


1 1 f\f\f\ 

11.900 


17.12 


O TOO 

3.722 


5.9 




1.975 


PM G016. 9-02.0 






3.508 




1.477 


7.537 


PM G018.6-02.2 


2.796 


17.56 




47.0 


1.802 


6.925 


PIN G019. 2-02.2 


3.370 


16.94 


O O f O 

3.852 


A C\ 

4.0 


1.778 


8.010 


PIN G019.4-05.3 


O A O O 

3.433 


1 C Of 

16.86 


A 1 1 A 

4.114 




on r 

0.865 


4.113 


PN G019.7+03.2 


3.208 


16.78 






1.414 


7.956 


PN G019.8+05.6 






3.576 


OO O 

28.0 


1 A Af\ 

1.449 


r oto 

5.872 


PN G023.8-01.7 


4.833 


17.39 










PN G023.9+01.2 






3.795 








PN G025.3-04.6 










1.630 


5.944 


PN G026.0-01.8 


7.225 


17.42 






1.912 


6.934 


PN G027.6-09.6 


3.825 


17.13 




0.6 


1.709 


1.660 


PN G031. 0+04.1 


4.238 


16.55 






1.575 


0.337 


PN G032.5-03.2 












0.175 


PN G032.9-02.8 


7.109 


17.05 


3.619 




1.455 


8.572 


PN G034.0+02.2 


4.755 


17.37 


3.077 




1.328 


6.880 


PN G038.4-03.3 






3.162 








PN G038.7-03.3 






3.722 


11.4 


1.894 


1.993 


PN G041.8+04.4 












0.135 


PN G042. 0+05.4 


21.790 


17.31 


4.046 








PN G042.9-06.9 


13.790 


17.05 


3.428 


0.3 


1.531 


5.184 


PN G044.1+05.8 












1.953 


PN G045.9-01.9 


11.190 


17.13 






-0.391 


0.977 


PN G048.1+01.1 


6.010 


16.66 


3.593 




1.309 


9.328 


PN G048.5+04.2 


13.460 


17.51 




66.0 




7.978 
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Table 4 — Continued 



Name 


Rr 


log Rpinj 


loe; N P 


1 A«Q« 
-"-4000 


Ne ratio 


EC 




fkpcl 


[cml 


fem - ^1 
i j 


[1/3=100] 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


PN G051.0+02.8 












6.412 


PN G052.9+02.7 


8.846 


17.72 






1.535 


6.759 


PN G052.9-02.7 


7.962 


17.06 








8.532 


PN G055.1-01.8 


20.580 


17.70 






1.879 


7.371 


PN G055.5-00.5 


6.597 


17.09 


4.114 




1.252 


2.061 


PN G058.9+01.3 


8.342 


17.43 






1.167 


2.389 


PN G059.4+02.3 


10.420 


17.37 






2.527 


6.327 


PN G059.9+02.0 


22.410 


17.11 










PN G060.5+01.8 


9.297 


17.22 






-0.055 


1.660 


PN G063.8-03.3 


19.380 


17.01 










PN G067.9-00.2 


7.743 


16.50 


3.787 




0.804 


1.526 


PN G068.7+01.9 


10.950 


17.41 






1.898 


5.166 


PN G068. 7+14.8 








2.4 




2.682 


PN G069.2+02.8 


13.670 


16.33 








0.153 


PN G077.7+03.1 






3.473 




1.044 


12.150 


PN G079.9+06.4 








77.0 




8.587 


PN G082.5+11.3 


12.530 


16.82 






0.895 


1.102 


PN G088.7+04.6 


11.620 


17.45 






1.846 


2.290 


PN G095.2+00.7 


10.110 


17.03 






1.158 


5.648 


PN G097.6-02.4 


14.730 


17.57 




7.0 


2.041 


5.935 


PN G104.1+01.0 


12.060 


16.88 






0.913 


6.476 


PN G107.4-00.6 












0.410 


PN G107.4-02.6 


14.200 


17.62 




92.0 




9.418 


PN G112. 5-00.1 


21.100 


17.68 


2.033 




0.833 


1.418 


PN G184.0-02.1 


15.270 


17.09 






0.012 




PN G205.8-26.7 












4.045 


PN G235.3-03.9 


15.260 


17.10 


4.689 






0.072 


PN G263.0-05.5 


12.990 


17.31 


3.762 


9.2 


1.915 


5.598 


PN G264.4-12.7 


13.390 


17.33 


3.949 




0.750 


2.723 


PN G274.1+02.5 












1.620 



Table 4 — Continued 



Name Rg log Rpn log N e I4686 Ne ratio EC 

[kpc] [cm] [cm" 3 ] [1/3=100] 

(1) (2) (3) " (4) (5) (6) (7) 



PN G275.3-04.7 


15.170 


17.33 




28.0 




5.872 


PN C278.6-06.7 


12.160 


17.23 


3.619 


6.0 


1.887 


6.354 


PM G285.4+01.5 


8.625 


17.27 


A 1 1 A 

4.114 




0.601 


5.656 


PM G285.4+02.2 






2.909 


78.0 




8.642 


PN G286.0-06.5 


16.060 


17.01 


O O A O 

3.348 




1 on 1 

1.391 


A A 1 r\ 

4.419 


PN G289. 8+07.7 


12.660 


17.39 




A r r\ 

45.0 




6.814 


PM G294.9-04.3 


10.210 


17.11 






-0.994 


0.661 


PIN G295.3-09.3 










1.362 


4.338 


PN C296. 3-03.0 


8.029 


17.12 


3.852 


19.0 


1.475 


5.374 


PIN G297.4+03.7 






3.299 




-1.256 


0.162 


PN G300. 7-02.0 


7.140 


17.21 


4.368 




1.002 


4.815 


PN G307.5-04.9 


6.532 


17.50 


3.787 


0.5 


0.034 


1.323 


PN G309.0+00.8 










0.773 


4.995 


PN C309.5-02.9 






3.302 




1.563 


6.934 


PN G311. 1+03.4 










0.289 


0.608 


PN G321.3+02.8 


5.232 


17.03 


4.647 




0.226 


3.213 


PN G324.2+02.5 


6.178 


17.37 


4.053 








PN G324.8-01.1 


5.434 


16.98 






1.352 


5.769 


PN G325.0+03.2 


5.935 


17.09 


3.366 




1.543 


6.530 


PN G325.8-12.8 


5.207 


17.49 






-0.057 


0.486 


PN G326.0-06.5 


4.562 


17.45 










PN G327.1-01.8 


4.366 


17.29 


4.114 




-1.617 


0.198 


PN G327.8-06.1 


13.910 


17.48 


3.677 




0.553 


2.362 


PN G327.9-04.3 






3.308 






5.062 


PN G329.4-02.7 


7.433 


17.47 


2.782 




1.702 


4.270 


PN G331.0-02.7 


4.750 


17.34 


4.339 




-0.703 


0.666 


PN G334.8-07.4 










-0.794 


0.688 


PN G336.3-05.6 


4.326 


17.36 


3.722 


56.0 


1.120 


7.424 


PN G336.9+08.3 












4.415 


PN G340.9-04.6 










1.902 


6.322 



Table 4 — Continued 



Name Rg log Rpn log N e I4686 Ne ratio EC 

[kpc] [cm] [cm" 3 ] [1/3=100] 

(1) (2) (3) " (4) (5) (6) (7) 



PN G341.5-09.1 


2.659 


17.50 


3.206 




1.889 


4.689 


Ti~l\T /"I O A O /( ill n 

PJN G343.4+11.9 


8.748 


17.51 




19.5 




5.401 


PJN G344.2+04.7 






3.619 




-0.020 


1.809 


PJN G344.4+02.8 












^7 /''no 

7.628 


ni\T f \ 1 ,1 no a 

PJN G344.8+03.4 






on/" 1 

3.306 




-0.695 


n A OT 

0.437 


PIN G345.0+04.3 














PIN G348. 4-04.1 






roo 

3.522 




1.427 


b.867 


T ") \ T f \ O Jn O f~\f~\ C~\ 

PIN G348. 8-09.0 










n ^ro 

0.753 




PJN G350.8-02.4 






A OTA 

4.250 




n on 

0.380 


1.737 


PJN G351.3+07.6 










1 nn a 

1.004 


2.871 


T~ii\T non m 

PN G351.9-01.9 






O A At\ 

3.449 




1 oon 

1.220 


5.665 


PN G352. 6+03.0 


1.464 


17.36 


3.779 




1.098 


4.838 


Tvt /"i or ^ 1 n ^ 

PJN G354.2+04.3 


2.556 


17.49 


3.533 




0.142 


0.661 


Ti"i\T /"i or ^ n 1 no r* 

PN G354.9+03.5 






3.141 




-0.027 


1.120 


PN G355.2-02.5 


3.563 


17.53 


3.895 




0.963 


5.252 


PN G355.7-03.0 


4.252 


17.41 


3.531 




1.351 


4.671 


PN G355.9+02.7 


4.518 


17.27 


4.101 




1.257 


6.313 


PN G356. 2-04.4 


0.540 


17.31 




4.5 


1.677 


2.439 


PN G356.5+01.5 


2.253 


17.56 






0.433 


1.854 


PN G356.5-03.6 






3.408 




1.173 


3.672 


PN G356.8+03.3 


11.350 


17.42 






-1.377 


0.211 


PN G357.1+01.9 






2.413 




1.324 


2.259 


PN G357.1-06.1 






2.879 


74.0 


0.998 


8.421 


PN G357.2+02.0 






3.762 




1.536 


9.364 


PN G357.6+01.7 


1.896 


17.32 


4.060 




1.390 


5.787 


PN G357.6+02.6 


5.437 


17.68 






1.359 


5.126 


PN G358.2+03.6 


1.449 


17.35 


3.583 


14.5 


1.675 


5.124 


PN G358.3+01.2 






3.839 




1.058 


11.772 


PN G358.5+02.9 












8.496 


PN G358.5-04.2 


0.237 


16.86 


4.114 




0.167 


2.430 
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Table 4 — Continued 



Name 



(1) 



Rg 
[kpc] 

(2) 



log Rpn 
[cm] 
(3) 



log 



cm 



(4) 



I4686 
[1/3=100] 

(5) 



Ne ratio 



(6) 



EC 

(7) 



PN G358.6+01.8 
PN G358.7+05.2 
PN G358.7-02.7 
PN G358.9+03.4 
PN G358.9-03.7 
PN G359.3+03.6 
PN G359.4+02.3 



3.723 
4.051 



17.36 
17.10 



3.792 
3.073 
3.058 

4.114 



1.605 



12.150 
6.066 



0.089 2.020 
-0.375 0.445 
1.337 7.758 



Table 5. Morphological distribution among dust types 



dust type 


N 


R 


E 


Symm. 


BC 


B 


P 


Asymm. 






% 


% 


% 


% 


% 


% 


% 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F 


5 




80 


80 


20 






20 


CRD 


24 


13 


61 


74 


4.2 


17 


4.2 


26 


ORD 


17 


12 


35.3 


47 




35.3 


18 


53 


MCD 


12 




50 


50 




25 


25 


50 



Table 6. Galactic bulge, disk, and Magellanic Cloud PN dust class distribution 



dust type Bulge, % Galactic, % Magellanic Clouds, % 
(1) (2) (3) (4) 



F 

CRD 
ORD 
MCD 



4 
11 

38 
46 



20 

28 
28 
24 



41 
52 
7 

N.A. 



